Method for Solving the Many-Body Bound 
State Nuclear Problem 



M. Fabre de la Ripelle, S. A. Sofianos*, and R. M. Adam 1 

Physics Department, University of South Africa, P. 0. Box 392, Pretoria 0003, 

South Africa. 

E-mail:* sofiasa@science.unisa.ac.za 



Abstract 

We present a method based on hyper spherical harmonics to solve the nuclear many- 
body problem. It is an extension of accurate methods used for studying few-body 
systems to many bodies and is based on the assumption that nucleons in nuclei 
interact mainly via pairwise forces. This leads to a two-variable integro-differential 
equation which is easy to solve. Unlike methods that utilize effective interactions, 
the present one employs directly nucleon-nucleon potentials and therefore nuclear 
correlations are included in an unambiguous way. Three body forces can also be 
included in the formalism. Details on how to obtain the various ingredients entering 
into the equation for the A-body system are given. Employing our formalism we 
calculated the binding energies for closed and open shell nuclei with central forces 
where the bound states are defined by a single hyperspherical harmonic. The results 
found are in agreement with those obtained by other methods. 

Key words: Nuclear structure, Hyperspherical Harmonics, Integro-differential 
equations. 
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1 Introduction 



The discovery by Jensen and Goeppert-Mayer of the shell model structure of 
nuclei [1] , was a great step forward in understanding the organization of nucle- 
ons in nuclei. Since the magic numbers occurring in the shell model can be ex- 
plained in terms of harmonic oscillator quantum numbers, with the occurrence 
of a large spin-orbit force, it was taken for granted that the nucleons move 
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inside an average one-body potential. Indeed, Goeppert-Mayer and Jensen [1] 
assumed that " • • • each nucleon moves in an average field of force V(r), of 
spherical symmetry, and independent of the exact instantaneous position of 
all the other nucleons. " In order to explain the existence of magic numbers, 
this one body field was completed by the introduction of a strong spin-orbit 
force. Some interesting remarks, however, have been put forward by de Shalit 
and Talmi about the validity of this interpretation [2], namely, that if such 
a self-consistent central field is found, then one should consider the residual 
interaction, which is not accounted for, as a perturbation. 

The shell model was introduced at the beginning of the fifties at the time 
where the nucleon-nucleon potential was not yet well known because the colli- 
sion energy was not large enough to detect the strong repulsive core needed to 
reproduce the S'-phase shifts beyond 300 MeV. This strong repulsive core gen- 
erates important two-body correlations in nuclei in such a way that Hartree- 
Fock methods applied with actual realistic nucleon-nucleon potentials produce 
a very small binding energy. Thus, the residual interaction responsible for the 
largest part of the binding energy cannot be considered as a perturbation. 
Of course if one considers that the free nucleon-nucleon potential is not the 
true interaction between nucleons in nuclei, the previous remarks do not apply 
but then one has to state it explicitly and indicate what interaction should 
be used instead. For instance in the self-consistent mean-field models [3] a 
Hartree-Fock method is used with Skyrme interactions. 

The shell model was first established for atoms but as pointed out by de 
Shalit and Talmi, there is an important difference between nuclei and atoms, 
namely, in atoms electrons have as a natural reference point the nucleus and 

. . .no such reference point exists in the nucleus, and at a first look it seems 
hard to understand why a central potential may form a good starting point for 
a nuclear model" [2]. 

Fifty years after the introduction of the shell model, the situation was still the 
same. After a lengthy discussion about the validity of the shell model, Talmi 
[4] pointed out that ". . . Today, as in 1949, the best proof for the validity of 
the shell model is the good agreement of predictions with experiments". The 
only well established property is that the quantum numbers associated with 
the magic numbers are those of the harmonic oscillator in a ground state 
corrected by a strong spin-orbit force. 

In order to find an alternative explanation, to that given by the Indepen- 
dent Particle Model (IPM), to interpret the magic numbers, let us start 
from the harmonic oscillator (HO) model. For an A identical-particle system, 
the HO potential, neglecting the parameter defining the strength, is propor- 
tional to the sum of the squares of the linear coordinates of the particles, 
Vffo( r ) ~ Hi x 1- This sum, when expressed in the center of mass frame in po- 
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lar coordinates in the D = 3(A — 1) dimensional space spanned by the particle 
coordinates Xi (i — 1, • • • , A), is nothing else but the radial coordinate, called 
hyperradius. The other space variables are the angular coordinates f2. Let f2 
and r be the polar coordinates fixing the position of the particles. For any po- 
tential V(r) which depends only on the radial coordinate r, such as for the HO 
potential where V(r) ~ r 2 , the rotational and vibrational motion described 
by a function of the radial coordinate are independent. The rotational motion 
in the physical three-dimensional space is described by spherical harmonics 
which are harmonic polynomials (HP) when the radial coordinate is equal to 
unity. In the .D-dimensional space, the rotations are similarly described by hy- 
perspherical harmonics which are HPs when the hyperradius is equal to one. 
These harmonics are independent of the shape of the potential V(r) since the 
rotational and vibrational motions are independent. Nevertheless, the speed of 
the rotation enters into the radial equation through the repulsive centrifugal 
barrier C(C + l)/r 2 , C = L + (D - 3)/2, where L is the degree of the HP 
defining the state in the D-dimensional space. 

Whatever the potential V(r) is, the ground state is obtained when the re- 
pulsive centrifugal barrier, i. e, the degree of the HP defining the state, is the 
smallest. For bosons where all particles can be in the same S'-state, the de- 
gree is L m = 0. For identical fermions where two particles cannot be in the 
same state, the wave function, which must be antisymmetric in the exchange 
of two particles, is the product of an antisymmetric HP and a function of the 
hyperradius r. This HP is a Slater determinant constructed from the indi- 
vidual polynomials sP^a^'^'Y^ '(u^) where s and t are the spin and isospin 
states and i refers to the rows and j to the columns for % and j running from 
1, • • • , A. The degree of this HP is the smallest when all the rrij quantum num- 
bers are used for the smallest possible values of 2nj + £j. It corresponds, as 
we shall show in the next section, to a HO Slater determinant in ground state 
and, therefore, to the filling of the states giving the minimum energy for the 
HO potential leading to the shell model. From this analysis it turns out that 
the " average field of force V(r) of spherical symmetrtf Goeppert-Mayer and 
Jensen speak about [1], which "represents the action of the other nucleons" is 
that part of the nuclear potential, called hypercentral, which is invariant by 
rotation in the D-space and thus a function of the hyperradius only. 

For the spin-orbit force one uses the property that a sum over all pairs of 
pairwise spin-orbit operators is equivalent to a sum over all particles of one- 
body spin-orbit operators [2]. When the spin-orbit interaction is the product 
of a hypercentral potential and a sum of pairwise spin-orbit operators, then 
this is equivalent to the product of a hypercentral potential and a sum of one- 
body spin-orbit operators [1] , but this time the potential is not the one of an 
independent particle model but a function of the collective radial coordinate, 
i. e, the hyperradius. The solution obtained from the wave equation where the 
potential is purely hypercentral does not contain any correlations since r is a 
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collective variable that describes only monopolar excitations when the radial 
wave function has nodes. The wave function has therefore to be improved 
when we have to deal with two-body potentials generating correlations. 

In this paper a method for introducing the correlations generated by the nu- 
clear potential is proposed. It leads to a system of coupled integro-differential 
equations taking the various operators of the potential into account. It is an 
extension to many bodies of accurate methods used for studying the few-body 
systems. 

In Sect. II we briefly describe the antisymmetric harmonic polynomials and in 
Sect. Ill the handling of the nuclear problem that neglects two-body correla- 
tions namely via the hypercentral approximation and the spin-orbit force. In 
Sect. IV we discuss the Integro-differential equation formalism while in Sect. 
V we describe in details the construction of the so-called weight and pseudo- 
weight functions associated with exchange operators followed by a description 
of the projection function in Sect. VI. The use of nuclear potentials is de- 
scribed in Sect. VII while a brief discussion concerning spin-isospin exchange 
generated elements is given in Sect. VIII. Our results are given in Sect. IX 
and our conclusions in Sect. X. Finally, details on the coordinates, and on the 
various coefficients needed in the formalism are given in appendices A, and B. 



2 Antisymmetric Harmonic Polynomials 

Any analytical function, including the wave function of a many-body system, 
can be expanded in terms of HPs or equivalently in terms of Hyperspherical 
Harmonics (HHs) multiplied by coefficients depending only on the hyperra- 
dius. Any such expansion starts from a HP of minimal degree L m called ground 
polynomial. When we have to deal with a hypercentral potential, i.e, with a 
potential invariant by rotation in the D-dimensional space spanned by the 
coordinates of the particles of an A-bodj system in the center of mass frame, 
this polynomial describes a ground state when the repulsive centrifugal bar- 
rier is minimal. When the potential is not hypercentral, the product of the 
potential with a ground polynomial generates new harmonic polynomials of 
higher degree. 

For the ground state and low excited states we expect, as was found for 
medium light nuclei, that the ground polynomial brings the largest contri- 
bution to the HH expansion of the wave function of the state under investi- 
gation. Since the other polynomials are obtained by starting from the ground 
polynomial, we define the state by the quantum numbers defining the ground 
polynomial. Ground polynomials for identical fermions must be antisymmetric 
in any exchange of two particles. In order to construct a ground polynomial 
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suitable for describing a system of A fermions, we introduce homogeneous 
polynomials 

<f) i (n,£,m)=x 2n+e Y ern (u i ) (1) 

of degree 2n + £ in terms of the polar coordinates (xi,u>i) of the ith particle 
where i — 1, • • • , A. 

In nuclear physics each nucleon i has a spin s and an isospin t and thus to each 
(f>i(n,£,m) polynomial a spin-isospin state st is associated. Let us construct a 
Slater determinant with the individual states ^i{n^£^m^)s\t\ where i defines 
the rows and j the columns for % and j running from 1 to A 

D=\\ 44Hnj^mj) || (2) 



If for each independent state (s, t, £, m) the determinant is constructed by using 
all the n from n = to n max (s, t, £, m), where n max is defined independently 
for each s, t, £, m, then this determinant is a harmonic polynomial. Indeed, if 
the Laplace operator 

A = X> 

i=i 

is applied to any column it decrease by 2 the degree of the homogeneous 
polynomial D without changing the quantum numbers (s, t, £, m). Since D has 
been constructed in such a way as to generate the lowest degree polynomial for 
the selected (s,t,£,m) quantum numbers, such a polynomial does not exist. 
Thus AD = and the antisymmetric homogeneous polynomial D is harmonic. 

Substituting the Laguerre polynomial L e + l l 2 (xl) for x 2n in fa one gets a new 
determinant 




If the same procedure of selecting the quantum numbers is applied, then ac- 
cording to the rule in which a determinant with proportional columns disap- 
pears, the new determinant, except for a normalization constant, is identical 
to the previous one. This new determinant multiplied by an exponential 

£>exp(-X>-/2) 
i 

is a harmonic oscillator Slater determinant [5]. This property will enable us 
to calculate two-body matrix elements by using HO Slater determinants. 

Instead of using spin-isospin and HO quantum numbers n, £, mi one can as 
well combine the spin s and angular momentum £ to generate a total angular 
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momentum j = £ ± ^ denned for each particle. Then the state of each nucleon 
is defined by the isospin t and the quantum numbers n,£,j,rrij where rrij is 
the projection of the total angular momentum j. This description will be used 
for open shell nuclei where the spin-orbit interaction contributes significantly 
to the binding energy. 

The removal of excitations of the center of mass motion occurring in the 
product of individual wave functions, is a difficult problem in the case of an 
IPM. As discussed by Talmi ([4], p. 51) "• • • it is only in the case of the HO 
potential that eigenstates can be written as product of an intrinsic function and 
a function of center of mass coordinates ■ ■ ■. Removing the contribution of the 
center of mass motion can be done only approximately • • • . It is still important 
to realize that such corrections are the price we pay for the great convenience 
of using shell model wave functions of nucleons moving independently in a 
potential well. " In the hyperspherical model the expansion of wave function 
starts from a ground HP translationally invariant. Indeed, by substitution 
of Xi — X (X for the center of mass) for X{ in Eq. (1) and according to 
the properties of determinants, Eq. (2), constructed according to the rule 

— * 

explained above, the X dependence disappears. Since the next terms to be 
included in the wave function are correlations described by functions of the 
relative coordinates r^ = Xj — Xj the overall wave function is translationally 
invariant and the center of mass problem is solved. 

Ground states are defined by ground HP where all available nucleons eigen- 
states of the shells are occupied up to the last one which, except for closed 
shell nuclei, are not fully filled. Excited states where one or several nucleons 
of the last shell are raised to the next shell are still described by a ground HP. 
But deep hole states where a nucleon is extracted from the core and raised to 
a higher shell, does not generate a ground HP and therefore such states can- 
not be described by our scheme. This results agree with the remark of Talmi 
([4], p. 7) that "• • • simple shell model states, like single nucleon or single hole 
states, if they lie at sufficient high energy, are not pure single hole states of 
deep-lying orbits have high excitations and may no longer be accurately de- 
scribed in terms of a single hole wave functions, the deeper the state of the 
missing nucleon, the higher the excitation energy of the hole state. The frag- 
mentation of single nucleon or single hole states at higher excitations impose 
practical limitations to the simple shell model picture. " 



3 The Hypercentral Approximation 

Let us start the description of our formalism by considering the case where 
correlations are ignored in the study of the A-body system, i.e, by discussing 
the hypercentral model first. The HO potential, neglecting the HO parameter, 
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V HO {x) 



A 

E 

i=i 



A 



X) 2 + AX 2 



(3) 



where X is the center of mass. The last term in the right-hand-side is irrelevant 
for isolated systems and therefore can be omitted. The other term is the sum 
of the square of the radial coordinates of the D — 3(A — 1) dimensional space 
spanned by the coordinates Xi, i = 1, • • • , A, of the A-particles in the center 
of mass frame expressed in polar coordinates (r, Q) . This sum defines the 
hyperradius r, 

r 2 = 2j2(x t -Xf = jj:(x l -x J ) 2 , (4) 

i=l i,j=i 



where the factor 2 is introduced so that for two particles r 2 = {x-y — X2) 2 . 

The hyperradius is invariant by rotation in the D- dimensional space. Here 
'hyper' refers to a D-dimensional space with D > 3. Since the HO potential is 
invariant by rotation in the D-space one may write the Schrodinger equation 
in polar coordinates (r, Q) where Q is the set of D — 1 angular coordinates on 
the unit hypersphere (r = 1) in the D-space which together with r fixes the 
coordinates x(xi, ■ ■ ■ ,xa) of the particles. Let V(r) be a potential invariant 
by rotation in the D-space. The Schrodinger equation, 

|-^^A i + V(r)- J B|* = 0, (5) 



in polar coordinate becomes 



m 



1 d D _, d L\n) 

,D-1 " " ^ 



dr 



dr 



+ V(r) - E \V = 



(6) 



where L 2 (Q) is the square of the grand orbital operator which is a gener- 
alization of the orbital momentum £(u>) in the 3-dimensional space to the 
D-dimensional space. It should be noted that in Eq. (6), the vibrational and 
rotational motions are independent. 

Let Hl{x) be a HP of degree L. Since it is a homogeneous polynomial, the 
r-dependence can be factorized out, Hl(x) = r L HL(Q). The Hl(Q) is a HH, 
i.e., the value of H L (x) on the unit hypersphere r — 1. From the assumption 
that H L (x) is a HP, with the definition of A in terms of r and f2 given in Eq. 
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(6), one obtains 



[L 2 {Q) + L{L + D- 2)]H L (n) = . 



(7) 



The wave function \1> can be written as a product of a HH and a hyper-radial 
function, 



1 



(8) 



where the hyper-radial function ui(r) is a solution of 

+ V(r) -E)u L (r) =0 



h 2 \ d 2 c(c + i) 



m 



dr 2 



(9) 



with C = L + (D — 3)/2 and where [L] stands for the set of quantum numbers 
including spin and isospin defining the state of grand orbital L. 

Whatever the potential V(r) is, the ground state is obtained when the repul- 
sive central barrier is the smallest. It corresponds to a HP of minimum degree 
L m . For a system of bosons where all particles can be in the S'-state this degree 
is L m = 0. For fermions, L is the smallest when HL m (x) is a Slater determinant 
constructed from the individual polynomials x 2n3+£j Y™ 3 (u>i), where i refers to 
the rows and j to the columns, when all the rrij are used for the smallest possi- 
ble values of 2rij +£j. As we explained in the previous section, this corresponds 
to a Slater determinant for a HO in a ground state. The minimum degree is 
the sum over the degree of the individual polynomials in i?[L m ](x) which is, 
as we have seen in the previous section, the sum over the radial and orbital 
quantum numbers of all the HO occupied states, 



3=1 



(10) 



The L m is related only to the rotation in the D-space and is independent 
of the shape of the hypercentral potential V(r). It is thus clear that in the 
hyperspherical scheme the "average field of force V(r) of spherical symmetry" 
Mayer and Jensen speak about, which "represents the action of the other 
nucleons" , is the hypercentral potential generated by the two body nuclear 
potential V(fij). This average potential is given by 

v(r) = A(A ~ 1] J H( Lm] (n)v(f tl )H [Lm] (n) dn (11) 
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where fij = Xi—Xj and V(fij) contains the spin and isospin exchange operators 
and in general all the components of the nuclear interaction V{fij) and where 

/ \H [Lm] (Q)\ 2 dQ = 1 . (12) 



It depends on the collective variable r and does not generate any correlations 
since when r is fixed, a variation in the distance between two particles does 
not affect the hypercentral potential and the corresponding wave function. 

To investigate the effect of a spin-orbit force we recall that it operates only on 
the last open shells. Following the derivation of de Shalit and Talmi we write 
for the spin-orbit operator for the pair {&% + Bj) ■ £%j where 

1 1 - - 1 

k = 2&~ f i) x - Pi) = + ~ 2 ^ j xPi + Xi* Pi) 

After elimination of the center of mass, the sum over all pairs of the spin-orbit 
operator gives 

A 

£ {Bi + Bj) • 4 = A^l-Si + S -L (13) 

i,j>i i=l 



where 

A A 

si = Bi/2, S = J2$ii L = J2^- 

i=l i=l 

The first term in the right hand side of Eq. (13) is the sum of one body 
spin-orbit operator occurring in the standard shell model and the last one, a 
collective spin-orbit operator, is a small perturbation. 

In the jj-coupling scheme, the spin-orbit operator does not modify the struc- 
ture of the Slater determinant defining the state and thus one has 

£ + a,) ■ 4 H [L] (Q) = C LS H [L] (Q) . 

i,j>i 

The coefficient Cls is determined by the filling of the shell of the determinant 
in the jj-coupling. 

The hypercentral potential, including the spin-orbit force Vg s (rij) (<7j+ Bj)- 
operating on a Slater determinant in the jj-scheme, is V(r) + C lsVls{t) ■ 
The ground state results from a competition between the repulsive centrifugal 
barrier h 2 /mC(C + l)/r 2 and the potential including the effect of the spin- 
orbit force through the coefficient C LS fixed by the filling of the individual 
states in the last shell. 
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In summary the quantum numbers defining the ground state are those of the 
HO in ground state for which the repulsive central barrier is the smallest 
corrected by the effect of the spin-orbit force. It is the model proposed by 
Goeppert-Mayer and Jensen but this time the potential is not a one body 
potential in which all nucleons move but a potential of the collective symmet- 
rical hyper-radial coordinate whose shape should be deduced from the realistic 
nucleon-nucleon potential. More information about the Hypercentral Approx- 
imation (HCA) where only the HH of minimal order L m is taken into account 
in the wave function can be found in Refs. [5-7]. We only note here that for a 
two-body potential V(rjj), the HCA provides a ground state binding energy 
for nuclei a little stronger but nearly identical to the one given by a variational 
calculation using a HO Slater determinant where the strength parameter is 
adjusted to give the lowest eigen-energy [5]. The HCA variationally provides 
an upper bound to the eigen-energy. A further improvement in the solution 
can be achieved with the inclusion of the two-body correlations generated by 
the two-body potential. 



4 Integro-differential Equations Approach 

For nuclear potentials with a strong repulsive core, the Hartree-Fock method, 
which reduces the interaction to a sum of individual potentials leading to 
an IPM, gives very poor results because the residual interaction responsible 
for the correlations is large and is not taken into account in the procedure. 
A Jastrow function might be added to take care of the correlations. In the 
method developed in this paper, one operates in the D-dimensional space 
spanned by the Jacobi coordinates and uses a radial coordinate system (r, Vt) 
where r is the hyperradius and Q stands for the angular coordinates at the 
surface of the unit hypersphere r — 1. 

When we have to deal with the one-body problem in polar coordinates and 
the potential in the physical three-dimensional space is not central, the wave 
function can be expanded in terms of spherical harmonics and the Schrodinger 
equation, projected on the spherical harmonics basis, generates a system of 
radial, coupled, second order, differential equations that have to be integrated 
to solve the problem. When we have to deal with many-body systems, treated 
in polar coordinates in the D-dimensional space spanned by the the coordi- 
nates of the particles, we have the same situation as for the one-body problem, 
except that in this case the space is larger. The interactions occurring in many 
body systems, for instance a sum of two-body potentials, are generally not in- 
variant by rotation in the D-dimensional space, and thus we have to deal in 
this space with a deformed potential. The expansion method can be applied 
in which HH must be substituted for Spherical Harmonics. Like for one par- 
ticle in a deformed potential, the Schrodinger equation is transformed into a 
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system of coupled, second order, differential equations in the hyperradius r. 
In this procedure we are facing a different problem namely that of degeneracy. 
While the (2£ + l)-degeneracy of the spherical harmonics for each orbital £ 
is moderate, leading to the use of a rather small number of significant terms 
in the expansion in the three-dimensional space, in contrast the degeneracy 
of the HH basis for a grand orbital L increases rapidly with the dimension 
D of the space leading to an intractable large number of significant coupled 
equations. 

A partial solution to this problem has been obtained by selecting the HH 
describing only two-body correlations, namely, the Potential Harmonics (PH) 
[6,8]. With this restricted basis, a good approximation can be achieved for the 
solution of the few-body problem. The rate of convergence can be improved 
by introducing functions limiting the number of significant coupled equations 
[9]. But for more particles, the number of coupled equations to be solved to 
obtain a good accuracy becomes again too large [10]. 

To overcome this difficulty the Integro-Differential Equation Approach (IDEA) 
has been proposed [11-13,7,14,15] in which the two-body correlations are taken 
into account. In the IDEA method, the Schrodinger equation is transformed 
into two- variable integro-differential equations, whatever the number of parti- 
cles is. The approximation is justified in saturated systems, like nuclei, by the 
rather low kinetic energy occurring between two particles because the volume 
occupied by the system increases proportionally to the number of particles. 
At this low kinetic energy, the pairs are mainly in S'-state and many-body cor- 
relations are expected to be rather small. It should be noted that the IDEA 
equation is identical to the Faddeev equation for three-particles when the in- 
teraction operates on pairs in an S'-state only. 

The capability of the IDEA to solve the Schrodinger equations for nuclei arises 
from the large component of the potential which is invariant by rotation in 
the D-space, i.e, from the hypercentral potential. The contribution of this 
potential to the binding energy balances the one of the kinetic energy in such 
a way that the binding energy in nuclei is mainly provided by correlations 
when a strong repulsive core (that give rise to large correlations) exists in the 
potential. In other words, we agree with Talmi [4] that • • the relevance to 
nuclear many-body theory is in realizing that two-nucleon effective interactions 
are all that is necessary to calculate nuclear energies. The two-nucleon matrix 
elements obtained by this procedure determines the structure of nuclei in their 
ground state and at low excitations. " In what follows we shall describe the 
IDEA formalism that takes into account two-body correlations and discuss 
various aspects of it. 
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4-1 Two-Body Correlations 



We assume that we are dealing with a system of A identical fermions of mass 
m, in particular nucleons. Our aim is to solve the Schrodinger equation 

(T + V(x)-E)V(a,x) = Q (14) 



where x = (x±, ■ ■ • ,xa) are the particle coordinates of an A-body system, a 
denotes the space independent degrees of freedom like spin and isospin, T is 
the kinetic energy operator, V(x) is the interaction potential, and E is the 
energy state of the system. The wave function can be expanded in terms of 
harmonic polynomials and the expansion begins with a ground polynomial 
H[L m ](x) of degree L rn for the state under consideration, 

oo oo 

¥(<*,*) = £ Y,H [L] (x)u [L] (r)= £ 5> w (fi)rV]M (15) 

L=L m [L] L=L m [L] 



where [L] is the set of quantum numbers, including the space independent 
degrees of freedom a, defining H[l](x). The notation H[l](Q) = Hm(x)/r L is 
used for the associated Hyperspherical Harmonics (HH) which is the value of 
H\n(x) on the unit hypersphere r — 1. Antisymmetric Harmonic polynomials 
are denoted by D^{x) (to remind us that they are constructed as Slater 
determinants) and the associated HH by D^{Q). 

To solve the Schrodinger Eq. (14) one needs to specify the properties of the 
interaction. When the interaction is a sum of one-body potentials V(x) = 
Yh=i V{xi) we have to deal with an Independent Particle Model (IPM) and 
the wave function is a product of individual states eigenfunctions of the one- 
body problem. When the potential is a many-body one, various methods have 
been proposed to obtain the solution of Eq. (14). In the most popular one, the 
state is described by the product of a Slater determinant constructed with in- 
dividual states adequately chosen and a Jastrow function which is the product 
of two-body functions Tii,j>i f( r ij) f° r describing two-body correlations origi- 
nating from a two-body potential. For most systems of identical particles like 
nuclei, the dominant part of the interaction is a sum of two-body potentials. 
The three-body potentials occurring next, can be reduced, except for a small 
contribution, to a sum of two-body potentials as well [16]. Thus in the present 
work we concentrate to the case where the interaction is a sum of two-body 
potentials only. 

In order to understand the structure of the solution, one starts from the HCA 
(discussed in Sect. 3) where only the ground harmonic, i.e, only the first term 
in the HH expansion of the wave function is taken into account. In this case, 
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the wave function, in the center of mass coordinate, is the product of a HH 
and a radial function 

*(a,2) = tf (r,n) = D [Lm] (n)- U5 ± m u (r) (16) 



where D[L m ](Q) is the ground antisymmetric HH for the investigated state. 

A projection of Eq. (14) on D[L m ](Q) for the wave function given by Eq. (16), 
results in the radial equation 



h 2 \ d 2 C(C + 1) 



m dr 2 r 2 



V [Lm] (r)-E 



«o(r) = (17) 



where C = L m + (D — 3)/2. The hypercentral potential Vj im ](r) is given by 
the integral 

v [Lm] (r) = J D* [Lm] (n)v(r,n)D [Lm] (n)dn (is) 



in terms of (r, fl) , the polar coordinates of x, taken over the surface of the unit 
hypersphere r = 1. In the three-dimensional space it would be the central part 
of the potential. When V(r, ft) is a sum of two-body potentials over all pairs 
and D\ Lm \{Vt) = H[ Lm ](Q) is normalized according to Eq. (12) then Eq. (11) 
holds. For the infinite hyperspherical well, the Coulomb, and the Harmonic 
Oscillator potentials the radial solutions Uo(r) are known analytically [5]. In 
the last case the Schrodinger equation can be solved either as a collective or 
as an independent particle model. 

To proceed beyond the HCA and obtain a solution of Eq. (14) with a better 
accuracy, one writes 

#(«,£) = = V (r,n) + *i(r,fi) (19) 



where *&i(r,Q) is the next improvement in the wave function which depends 
on the structure of the potential. Here we are interested for potentials which 
can be written as a a sum of two-body interactions 

V(x) = ]T V(r tJ ) , r l3 =x % - x 3 (20) 

i<j 

Using Eq. (19), Eq. (14) becomes 

(T + =E^i)(*o + *i) (21) 

i<j 
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Omission of the correction term on the right hand side, implies a structure of 
of the form 



* 1 = * ]TF(%r) (22) 

i<j 

where F( two body amplitude. Thus, for nuclei with central poten- 

tials we may write 

*(£)=£ [Lm] (fi)]TF( W ) (23) 

i<j 

in terms of the two-body amplitude F(rij,r) where the pair is assumed to be 
in an S'-state. 

In order to find the amplitude F(rij,r) we have two options: Either to solve 
the amplitude equation 

(T - E)D [Lm] (n)F(r ij: r) = -V(r tJ )D [Lm] (n) £ F{r kU r) (24) 

e,k<e 

or to solve the Schrodinger equation 

(T -£)*(£) = - £ V(r t ^(x) (25) 

j,i<3 

where now \P is given by Eq (23). The first option, Eq. (24), where the sum 
over all pairs reproduces the Schrodinger equation, leads to the IDEA [11- 
13,7,14,15], while the second leads to the Variational Integro-Differential Equa- 
tion (VIDE) [17,18]. It should be noted that for three-body the IDEA becomes 
a Faddeev-type equation while the VIDE does not have any counterpart. The 
solutions are not the same, as the VIDE is variational while the IDEA is not, 
but the latter is easier to solve. When the potential V(ry) operates only on 
pairs in 5-state the Faddeev equation provides the exact solution but when 
the potential is local and operates on all orbitals, the VIDE gives far more ac- 
curate solutions. In order to obtain the same accuracy, a few coupled Faddeev 
equations have to be solved. In both cases the wave function has the structure 
of a sum of two-body amplitudes. 

The correlations are described by a product when Jastrow functions are used. 
However, the product structure of the Jastrow function commonly used, does 
not have the correct structure of the solution for three or more particles with 
an S'-projected interaction. In contrast, in the IDEA formalism we only assume 
that a sum of two-body amplitudes solution of Eq. (24) is sufficient to provide 
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us with a good solution although many-body correlations are not taken into 
account. 



4-2 Reduction of the IDEA Equation 



We assume that the wave function is given by Eq. (23), i.e, by a product of 
an antisymmetric HH and a symmetrized two-body amplitudes. Our aim is 
to solve Eq. (24) to obtain the two-body amplitude F(rij,r) associated with 
the ground polynomial D[ Lm ](Q) defined by the occupied states in a Slater 
determinant. In order to have an equation for r\j only, for instance for a 

— * 

reference pair = £ N with N = A — 1, i = 1, and j = 2, we have to eliminate 
the part which in F(tm,t) depends on the other Jacobi variables i < N, 
which are related to many-body correlations (For the A-particle coordinates 
in the center of mass frame, see Appendix A). For this purpose, one projects 
F(rij,r) for 7^ (1,2) on the pair (1,2) in an S -stake. Let us call V c the 
operator projecting the connected pairs like (1, i) and (2, i), i — 3, • • • , A and 
Vd the disconnected pairs like for % and j > 2. Then the total projection 
operator is 



, , (A-2)(A-3) . 



V° = 2(A - 2)Vl + ^ § '-VI (26) 



where V® operates on one connected pair, e.g. the pair (2,3), and V® on one 
disconnected pair, e.g. the pair (3, 4). Since our aim is to calculate correlations, 
we isolate on the left hand side the hypercentral part of the potential and in 
the right hand side only that part of the potential generating correlations. 
Then Eq. (24) becomes 



(t+ A{A 2 % [Lm] (r) - E^D [Lm] (n)F(r 12 ,r) = -<V(r 12 ) - V [Lm] (r) 
xD [Lm] ({l)! [ F(r 12 ,r)+(A-2) [2V° c F(r 23 , r) + ^^°F(^4, r)] } (27) 



where the functions within the bracket { } depend only on r 12 and r. Eq. (27) 
is written for the pair = (1, 2). It is noted that by summing over all pairs 
one recovers the Schrodinger equation (25). 

The two variables equation in r and z we are looking for, which determine the 
two-body amplitudes, is obtained by multiplying Eq. (27) at left by Df L AQ) 
and by integrating over the surface element df2jv_i (see Appendix A). First 
we require the integral 

J D\ Lm] (Q)D [Lm] (fl) dlW = W<2, (z, U ) , (28) 
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taken over all angular coordinates fijv-i excluding z = cos 20 with r^/r = 
cos0 and u, the angular coordinate of r^, and thus this integral is a function 
of z and uj. We may define the weight function W[L m ](z, uj) for the state [L m ], 

W [I ^ ] (z,u>) = W (z)W™ ] (z,u>), (29) 
where 

W (z) = T D I\\ - zf D -^ 2 {\ + z) 1 ' 2 , (30) 

in such a way that for normalized D[ Lm ](Q) (see Eq. A. 8) 

+i 

J D* [Lm] (Q)D [Lm] (Q)dQ = J W [Lm] (z,u;)dzdcj = 1 (31) 



where dQ is the surface element on the unit hypersphere given by dQ = 
Wo(z) dzdudQw-i. The construction of W[L m ]{z,u) is discussed in the next 
section. We only note here that for bosons in ground state [L m ] = 0, 

D [0] {U) = Y [0 ] = 2?rD/2 , 

and W[o](z,u) reduces to Wo(z) except for a normalization constant. 

To proceed we must calculate the kinetic energy and the projection terms. 
For the former term we use the relation (A. 12) with the center of mass at rest 
omitted, 

h 2 h 2 A_1 
r = ~ A e = ~EVl, i = l,"-,A-l 



m m . 

i=i 



and set 



F ^s,r) = -^ m P(z,r). (32) 



Therefore, we have to calculate 



where 



1 d D _ x d L 2 (Q) 



~D-1 Q r Q r 



with the grand orbital operator being written as 
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W^[ ] (2) <9z 1 J 1 + 2 

+ derivatives in Qn~i coordinates 



(33) 



Thus, multiplying Eq. (27) from left by D* Lm ](fi), integrating over the surface 
element dfijv-i, and taking into account the relation 

2D iL m] m§- z D [Lm] (Q) = A [2? rLm] (n)£> [£m] (n) 



one gets the basic equation [11-13] 



h 2 



m 



d 2 | C{C + 1) 



dr 



+ ^^ 1 V [Lm] (r)-E ] jP(z,r) 

= - (V(ry/(1 + z)/2) - V [Lm] {r)) [p(z,r) +V°P(z,r) 



(34) 



where V°P(z, r) is the projection of the two-body amplitudes for all the pairs 
for 7^ (1)2) on the reference pair (1,2) in S'-state [13] (see forth 
Sect. 6). 

Eq. (34) can be easily solved directly as a two-variable integro-differential 
equation to obtain the binding energy of the nucleus under consideration. 
However, it is desirable to use also the adiabatic and the more accurate uncou- 
pled adiabatic approximations which can provide us not only reliable solutions 
but also the eigen-potentials from which further physical information, such as 
low energy scattering states, can be extracted. 



4-3 The Extreme Adiabatic Approximation 



In most systems the energy contained in the rotation in the D-dimensional 
space is very much larger than the one in the radial motion. Indeed, in nuclei 
the radial energy of the ground state is of the order of half the monopolar 
excitation energy (breathing mode), i.e, about lOMeV, compared to the sev- 
eral hundreds or thousands of MeV generated by the rotation which can be 
estimated from the Fermi gas model giving the average kinetic energy per par- 
ticle E kin /A = 28.7r^ 2 fmMeV, i.e, ~ 23.7MeV for r = 1.1 fin . This means 
that the rotation and the vibration are nearly decoupled and thus one can 
use an adiabatic approximation in which one freezes the r-motion and solve 
the rotational motion equation to obtain for each r an eigenpotential that 
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is subsequently used in the radial equation [19]. More specifically, one writes 
P(z,r) = P\{z,r)u\{r) and solves 



4 r 

= [V(r^/(l + ^)/2-V [Lm] (r)] 



P A (^,r)+P P A (z,r) 



(35) 



for fixed r to obtain the eigenpotential U\(r) which is then used in the radial 
equation 



d 2 £(£ + 1) 



dr 



+ ^4r^% m ]« + ^(r)-pU A (r)=0 (36) 



to obtain the total energy E and the radial function u\(r) [11]. The Eq. (35) 
and (36) constitute the so-called Extreme Adiabatic Approximation (EAA). 

The wave function, within this approximation, is obtained from P\(z,r) and 
u\(r), 



^eaa(x) 



1 



r 3A/2-2 



Mr) E Px(2rl/r 2 -l,r). 

i<j<A 



(37) 



When the hypercentral potential Vjx, m ](r) is ignored, in Eqs. (34-36) for A = 3, 
one finds the Faddeev equation for S-state projected potentials [20,21]. For 
A > 3 it will be called S-State Integro-Differential Equation (SIDE) and 
should be applied when we are dealing with S'-state projected potentials since 
V° in Eqs. (34-36) projects any other pair(i,j) on the pair (1,2) in S'-state. 
In the IDEA only the residual interaction generating correlations operates on 
pairs in S-state. 



4-4 The Uncoupled Adiabatic Approximation 



In order to separate the two- variable integrodifferential equation (34) into two, 
one-variable, equations (35 and 36) we assumed that the amplitude can be 
written as a product P(z,r) = P\(z,r)u\(r) and the EAA is obtained when 
the variation of P\(z,r) with respect to r is neglected. It provides a lower 
bound to the eigen-energy [22,21]. A further improvement in the binding can 
be achieved if such a variation is taken into account. 

The new equation is obtained by substitution of P\(z, r)u\(r) for u(r) in (36) 
and by taking the derivatives of P\(z,r) with respect to r into account. One 
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first normalizes the P\(z,r) according to 



+1 

(Px\Px) = J P 2 x (z,r)W [Lm] (z)dz = l 



(38) 



which leads to the derivatives 

dPx , . . _ , d 2 P : 



A 



dr 



0, (P; 



A 



dr 2 



A ) = _ ( dPx | dP A) 



dr dr 



(39) 



Multiplying Eq. (36) with P\(z,r)W[L m ](z) from right and integrating over z, 
we obtain the new equation for u x 



m 



d 2 £(£+1) 



dr 



^4(^4 1) , . 
+ J V [Lm] (r) 



+ U x (r) + -( 



h 2 , dP x . dP x 



m dr dr 



)-E\u x (r) = 0. 



(40) 



The solution thus obtained is called the Uncoupled Adiabatic Approximation 
(UAA). We note that the extra term introduced in (40) is always negative and 
the new effective potential 



(41) 



always provides an upper bound to the eigen-energy. 



4-5 A Variational Equation for the IDEA 



The IDEA Equation (34) is a two-variable integro-differential equation where 
the integral part comes from the projection V° of all pairs on the 

reference pair (1,2). The two-body amplitude P(z,r) is a solution of an 
integro-differential equation which must fulfill certain asymptotic conditions 
for z = ±1 and r = and r = oo. A variational equation can also be obtained 
in the adiabatic approximation by multiplying Eq. (35) by the weight func- 
tion, normalized to unity, and by integrating over z in the range [—1, +1]. The 
normalization is 

+i 

J P x (z,r)W lLm] (z)dz = l (42) 
-l 
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implying that the hypercentral part of P\(z,r) is taken to be 1. It is the 
dominant contribution in the HH expansion of P\(z,r). The kinetic energy 
term in the left hand side disappears and the equation for U\(r) becomes 



Ux(r) = J [V(ry/(l + z)/2-V [Lm] (r) 
-l 

P(z,r)+V°P(z,r)] W [Lm] (z)dz 



(43) 



Thus one may search for the solution P\(z,r) for which U\(r) is a minimum 
which is then introduced in the radial equation (36) to calculate the ground 
state energy. One can also control the quality of P x (z,r) solution of Eq. (35) 
by calculating the eigen-potential U\(r) from (43). This should be identical 
(within numerics) to the one obtained from Eq. (35). 



5 Construction of the Weight Function 



The HO wave function can be written either in the .D-space as a product of a 
HP D\L m T,(x) and an exponential in the hyperradius or as a Slater determinant 
constructed from individual HO eigenfunctions. The weight function can then 
be obtained by identification of the Fourier transform in the relative coordinate 
fij of |D[l](x)| 2 calculated first from the HP representation in the .D-space and 
then from the standard HO Slater determinant. 

Let D^(x) be the normalized Slater determinant describing the states of an 
A-nucleon system 



where i refers to the rows, j to the columns (i, j = 1, • • • , A) and where s and 
t denote the spin and isospin states. The normalized HO eigenfunction for one 
particle in polar coordinates (x, oj) is 

1/2 

Y e m (uj)(x/bYL e n +1 / 2 ((x/b) 2 ) e -(^ 2 / 2 (45) 



2n! 



6 3 r(n + £+f) 



where b is a parameter related to the size of the system. The determinant 
describing the HO ground state of the A-nucleon system is constructed in 
such a way that for a given set of quantum numbers s J ,P ,£j, m,j all the quan- 
tum numbers rij, for rij running from to a maximum value rij{s\ P, £j, rrij) 
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chosen independently for each set, are filled. According to the properties of 
determinants it can also be written as 

Df L f(x) = C L || s^xT^Y^^) || e-<£f*?>/«3 (46) 



where g 2 = 2b 2 , £f x 2 = r 2 /2 + AX 2 , X being the center of mass coordinate, 
and Cl is a normalization constant. This determinant is a translationally 
invariant ground HP of degree L = Y^f =1 (2rij + 

One may factorize r L and rewrite Eq. (46) as 

T 2r 2L e"^o 1 1/2 T 2e- R2 / fo2 1 V2 



where R 2 = AX 2 and Dm (Q) is a normalized HH of grand orbital L. Indeed, 

■)H0 



by integrating \Df[f(x)\ 2 over the 3A-dimensional space with 



d 3A a; = d 3 ^ • • • d 3 x A = dttr - 1 dr d 3 R (48) 



one finds 



/ \D?° (x)\ 2 d 3A x = J \D [L] (n)\ 2 dn = 1 . (49) 



The integral f \D[ L ](£l)\ 2 dil N _i (occurring in the weight function) can be 
obtained from the Fourier transform 

V(y,u k ) = (D?°W\e rk <»\Df°(x)), £ N = f 12 (50) 



calculated by two different methods. Since r L D^(Q) is a polynomial of degree 
L homogeneous in the Jacobi coordinates i — 1, ■ ■ ■ , N(= A — 1), then 
|r L .Drn(n)| 2 is an even homogeneous polynomial and the integral over Qn-i 
in Eq. (28) is a homogeneous polynomial in the £jv and p 2 = r 2 — £^ variables 

/ \r L D [L] (Q)\ 2 dn N ^ = £ ([L]\n,i}Y e (co)e N n+£ p 2Lm - (2n+e) (51) 

i even 



where ([L]\n,£) are coefficients that have to be defined. The volume element 
for the coordinate system (fl N _i,p, £ N ) in the D = 3N = 3(A— l)-dimensional 
space is 

d 3JV £ = dQ N ^ p D ~ 4 dp£ 2 N d£ N du (52) 
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where ui is for the angular coordinates of £at- Since the operator exp[ik ■ r*i 2 ] 

— * 

is independent of the center of mass, the integral over R in (50) gives 1 and 
the Fourier transform becomes 



£ even 



2(n+l)+^ 

dl^ 



/ \ 2L+£>-4-(2n+£) / \ 

*(*■) 4-) < 53 » 

\QoJ \QoJ 

One substitutes for the plane wave e ik ' x , (x = £ N ), the expansion in the polar 
coordinates (cu, x) of x, 

e &* = (2n)V 2 ^YrMYr(u)-}=Je +1 /2(kx) . (54) 



where (k, uj^) are the polar coordinates of k. After integration over the (p, £jv) 
variables, the Fourier transform becomes 

3/2 

1 + L> I Z ) n ,i 

x r(L + (D- 3)/2 - n - £/2)y<V*)^ +1/ V) e"" 2 (55) 
where y = kb/\/2 = kg /2. 

The Laguerre polynomials L"(x), x = y 2 , constitute a complete orthogonal 
polynomial basis associated with the weight function x a e~ x . Thus, multiplying 
Eq. (55) by Y® (cOk)y £ L^ 1 ^ 2 (y 2 ) and integrating over y in the range < y < oo 
and over Uk we obtain the expansion coefficient occurring in Eq. (51) 



OTn, €) = (-!) 



H2 ± T(L + D/2) 



tt 3 /2 F[L + (D- 3)/2 - {n + £/2)]T{n + £ + 3/2) 



x 



J V(y,u k )Y e °My e L^ 2 (y 2 )y 2 dydu k (56) 



According to Eq. (55), the Fourier transform is a polynomial in Y®(ujk)y 2,/+e 
which can also be written as 

V(y,u k ) = v^£(-l) A/2 ([L]M)nV^ +A e~ y2 (57) 
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where ([£]|^, A) are new coefficients to be determined. Introducing (57) in (56) 
and using the analytical expression 

J tfW Llt 1/2 {y 2 ) e^V dy = t^L h F(v + £ + 3/2) (58) 



we find the relation between the two kinds of coefficients ([L] \n, £) and ([^]|^, £) 



tt v ' T(L + (D -3)/2-n-£/2) 
x£^W + £ + 3/2)([L]|M)- (59) 

u>n VV 

Introducing (59) in (51) with £jv = rcos0, p = rsin0, and z = cos 20, 
dividing by r 2L , and using the definition of the Jacobi polynomials P" ,/3 (x) 
[23] one obtains 

2 (-l) n n! 



/ | D [L] (Q)\ 2 <m N -i = -r(L + D/2) 



r(L + (D-3-^)/2) 

I even 

x([L] \n, i)Y e °(uj) (cos0)'(sin 0) 2L - 2n ^P T r n ~' /2,/3 (cos 20) (60) 
where a = L + (D - 5)/2 and (3 = i + 1/2. 

Since i is even while sin 2 = (1 — z)/2 and cos 2 = (1 + z)/2 the sum over 
i and n generates a polynomial in z. We have now two expressions at our 
disposal for occurring in the weight function (29) according to whether 

one uses {[L m ] n,£) or ([L m ]\n,£) coefficients. In the first case 

W™(z,u>) = ± T,([Lm]\n,£) Y?{u){l + z)^ (1 - z )^-n-e/2 (61) 



and in the second 



^ (z u) - J_l V (-2)"n!r(L m + D/2) 

^"■1 U ' w) " 2^ tt ^ r(L m + (D - 3 - £)/2) ( [ ml |n ' £) 

X (1 - + ^y*^) pa-n-W+1/2^ (g2) 

In order to find explicit expressions for the coefficients occurring in either (61) 
or (62) we have to calculate the Fourier transform, Eq. (50), where D^(x) 
is the IPM representation of the Harmonic Polynomial D^(x), Eq. (44). For 
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this purpose we first write 



and then develop the normalized Slater determinant -Djx^j (^0 i n Eq. (44) with 
respect to the first two rows for xi, and x 2 



V ~ ^ i,i<j 



(63) 



dij(xi,x 2 ) 



where 

siiiV'i(^i), sUiV'j(^i) 

s* 2 ? 2 ipi(x 2 ), siti^j(x 2 ) 

Since -D^° and -D^>° do not contain the same HO individual state except for 
i — k and j = t we have 



(64) 



Integration over all Xj for i > 2 leads to 

/ \d*° (x)i 2 d 3 x 3 • • • d 3 x A = \d l0 {x u x 2 ) 



(65) 



where 



and 



Therefore, 



/ 



d lj (x 1 ,x 2 )- j.^ 
|<%(£i,£ 2 )| 2 d 3 ^ d 3 :r 2 = 2. 



The djj contains identical columns and thus it vanishes. Further, since d y - is 
invariant by exchange of % and j, the sum over % and j > % can be transformed 
into a sum over all states independently of their position in the determinant 



I d^fa, x 2 ) | 2 = i I ^2) | 2 



(66) 



For simplicity, we introduce a single set of quantum numbers ( for spin and 
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isospin states, labeled according to 



C — +| f° r proton spin up 

C — +| f° r proton spin down 

C — — | for neutron spin up 

( = — | for neutron spin down 

where we count 1 or -1 for proton and neutron respectively and |, — | for 
spin up or down and we add the values to obtain (. In order to calculate the 
Fourier transform we do not have anymore to refer to the columns but only 
to the spin-isospin and space occupied states. Then a sum over all pairs of 
occupied states is substituted for the sum over % and j > i in Eq. (65). 

Let us consider the pair of states ( a \a > and £p \/3 > where \a > and \/3 > 
stand for the HO eigenfunctions (45) with quantum numbers (n a , £ a ,m a ) and 
{np,tp,mp). Then one may distinguish two cases in the calculation of the 
Fourier transform 

F Q/3 = (^(12)1 e *<*-*> |^(12)) . (67) 



In the first case one defines the direct term 

D aP (k) = (a\ e il£ \a)(p\ e'^lP) + {k -> -k} (68) 

and the exchange term 

E a/3 (k) = (a\ e Ss \P){P\ e- iS "*\a) + {k ^ -k} (69) 



where k — > —k means the first term in which k is changed to — k. These Fourier 
transforms are real and even functions of k and thus they contain Y^{u) for 
even I only, where I > is related to the deformation of nuclei. 

In Nuclear interactions spin and isospin can be exchanged between two nucle- 
ons % and j. The exchange operators are Pff where e = without exchange 
while e — a, r, <jt where spin, isospin or spin and isospin are exchanged re- 
spectively. These operators are traditionally known as Wigner (e = 0, Bartlett 
(e = a), Heisenberg (P^ = —P[j, and Majorana (Pff = —P aT ), exchanged 
operators. In the above, the Fourier transform has been calculated without ex- 
changed operators i.e for e = only. This case corresponds to the Wigner force 
for which the weight function is given by Eq. (29) with (61) or (62). Therefore, 
we have still to calculate the so-called pseudo-weight functions associated with 
the operators Pfj, where e — a, r, or ar. 

To evaluate the contribution brought by the direct and exchange terms D a p 
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and E a p, to the Fourier transform 



F<$ = (^(12)| ^-^Pg\d a/3 (12)) 



(70) 



including the exchange operator P[ 2 for e = 0, a, r, err, we have to distinguish 
four cases: 



1) The spin-isospin states ( a and (p are the same, i.e, ( a — (p = 0. 

2) The spin-isospin states are different, i.e, ( a ^ (p but the P-^ operator does 



not modify the spin-isospin states. 

3) The operator P-p exchanges the two spin-isospin states. 

4) The operator pff generates new spin-isospin states orthogonal to the orig- 
inal ( a and (p states. 

In terms of the Kronecker symbol 5 a b = 1 for a = b and zero otherwise, the 
contribution in the first case is 



in the second 
in the third 



k^D^ - E a p), 



while in the last is zero. Applying these rules one finds that for the Wigner 
operator P^ = 1 we have the Fourier transform 



F a(3 = D a/3 - SCcCpEap. 



(71) 



Instead of considering the Pf 2 and P{ 2 operators separately, it is more conve- 
nient to use the combinations 1/2(P 1 CT 2 ± P[ 2 ) leading to the expression 



FS = \{l ~ % + c,)o){(l + k^)\(^ + 5 £T ) 



D 



a,f3 ~ E a ^ 



- (-1)^(1 -5^)^-5^) 



D at p + E a> p 



(72) 



where e = a or e = r. Finally, for the spin-isospin exchange operator or we 
have 



a/3 



kd D aP ~ E, 



a/3 ■ 



(73) 



From the above analysis we see that in order to calculate the weight and 
pseudo-weight functions associated with the operators occurring in realistic 
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nuclear potentials one needs to know the Fourier transform of the product of 
two individual HO eigenfunctions. Indeed, since we have to deal with a two- 
body weight function, which depends on z = 2rf-/r 2 — 1 (here = (1,2)), 
the Fourier transform of the product of two individual HO eigenfunctions 
occur twice for the states \a > and \/3 > in the calculation of the direct and 
exchange terms in Eqs. (68) and (69). 

5. 1 Weight and Pseudo- Weight Functions: General Case 

Let us define the Fourier transform 

/ r na/a , ma e^, W d 3 x = H e**|/J) (74) 

where \at > and \/3 > refer, respectively, to the HO eigenstates with quantum 
numbers n a ,£ a ,m a and np,£p,mp in the HO wave function (45). The Fourier 
transform is the sum over the quantum numbers n and t 

(a| e iU - s \(3) = y/^^2(i) t (a\n,£\P)Yp(uj k )y 2n+l e~ y2/2 (75) 

where m = mp — m a . The coefficients (a\n,£\P), normalized according to 
(a|0, 0\P) = S a/ 3, are given in the appendix. For j3 — a the simplified notation 
(a\n, £\a) = (a\n,£) is also used with (a|00) = 1. 

The calculation of the Fourier transform (57) depends on the direct and ex- 
change terms D a p(k) and E a p(k). However, the question of normalization must 
be addressed first. In Ref. [24] it was shown that for closed shell nuclei the 
coefficient ([-^]|^, A) occurring in (57) and in the weight function written as in 
Eq. (62) are simple and most of them integers and the normalization in Eq. 
(49) is A(A — l)/2 instead of I. To obtain afterwards the normalization to 1 
obviously one divides by A(A — l)/2. To obtain these coefficients one must 
calculate the Fourier transform for the expression 

F=\E( d ij\e iHSl - S2) \ d ij)- (76) 

i,3 

where the factor 4 stems from the normalization to A(A — l)/2 of this expres- 
sion for k = and from the fact that da = dij = —dji with (dij\dij) = 2. 
Using Eq. (75) and the expansion of the product of two spherical harmonics 

JTWK) = EPTOW) >TK)> ( 77 ) 

A 
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with n = mi +77^2, the direct term, Eq. (68), for the couple of states \a > and 
\(3 > with mi = m 2 = becomes 

D af} (k) = ^Y.(- l ) X ' 2d ^^\v ) \)Y^k)y 2v+X e~y 2 (78) 



where the direct coefficient 

d(a,/?|i/,A) = VSF ^(-^^"^("Ni^O^K^^lF/J^) (79) 



for v = ni+n2 + (^i+^2 _ A)/2 where £i, £2, and A are even and ^1 < A < 
£1 +£2- A factor of 2, stemming from the two terms in (68) for k and k — > —A;, 
has also been taken into account. One notices that for n x = n 2 = £\ = £2 = 

— * 

for which v = X = the normalization d(a(3\00) = 1 leads to D a p{k = 0) = 1. 
A similar procedure is applied to obtain the exchange term 

E aP {k) = £(-l) A / 2 e(«,/3k A)F A °K)^ +A e^ 2 (80) 



where the exchange coefficient 



e(a,(3\v,X) = V^ £ (-l)^^ 2 ^^!^, 41^(^1^, £ 2 |a) 

nl £j 

™2 ^2 

x (y A °|y-|y-> (81) 

for v = rii + 112 + (£1 + £2 — A)/2, A even, mi = mp — m a , m 2 = —mi and the 
normalization e(a,(3\0,0) = 5 a/ 3 

Let us come back to the Fourier transform, Eqs. (67) and 

f Vf L f*(x)e^-^Df L f(x) d SA x = - ]T F aP (82) 

all occupied 



for the normalization (76) where the wave function is normalized to A(A— 1)/2 
as explained above. Using (78) and (80) one finds, by identification with (57), 
the coefficient for e = i.e for = 1 in Eq. (71), 

{\L]\uAT ] = \ E {d(a,(3\u,X)-5 Ca ^e(a,(3\u,X)} (83) 

all occupied 
states 
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where the sum is taken independently over all HO occupied states in the Slater 
determinant D^(x) defining the state under consideration. This coefficient, 
introduced in Eq. (62), defines the weight function 



W lLm] (z,u>) = W^Jz,cv) = W {z)W^M 

according to Eq. (29). The upper index (0) is to recall that we are dealing 
with a Wigner potential. 

For the pseudo-weight function corresponding to the spin, isospin, and spin- 
isospin exchange operators P^ respectively for e = a, r, ot one uses Eq. (72) 
and (73) leading to the new spin and isospin coefficients 



all nrnminH ^ 



all occupied 
states 



+ ^ a ( P )^ea + S eT )^d(a, \) -e(a,P\v, A) 
_(_l)C+fc(i _ 5 Ca(: /-(5 ea - 5 eT ) [d(a, j3\v, A) + e(a, (3\u, A)] } 

for the Bartlett (e = a) and Heisenberg (e = r) forces and to 
1 



£ 



all occupied 
states 



x d(a, (3\v, A) — e(a, (3\u, A) 



(84) 



(85) 



for the Majorana force (e = ar). The weight and pseudo- weight functions are 
thus defined by 



W$Jz,u;) = W Q {z)W{fa{z,u) , e = 0,a,r,ar 



(86) 



where in Eq. (62) one uses ([L m ]|n, £) e instead of ([L m ]|n, €) and with normal- 
ization 



/ 



£>K(*)r<i*" 



/ 



(87) 



5.2 Weight Function for Coulomb Potential 



The Coulomb interaction is of Wigner type i. e, it does not contain any spin- 
isospin exchange operator, the interaction being considered only between charged 
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particles, in our case the protons. The projection operator applied to d a p{l, 2) 
is 

^^(1,2) = J(l + sgnCa)(l + sgnC/9K/3(l,2) (88) 
which cancels the neutron states. 

In the Fourier transform (70) with e = c (for Coulomb) only the first two cases 
contribute. Using the same procedure as for the other projection operators 
the two-body coefficients are given by 

m\^f = \ E (l + sgn(Ca))(l + sgn(^)) 

all occupied 
states 

x{d(a,/3|i/,A)-<J faf<J e(a,/3|i/,A)} (89) 

The hypercentral potential for the Coulomb interaction can be be obtained 
directly from the Fourier transform of the Coulomb weight function [24] with 

«(*) = ^ <*» 
the result being 

^ w = I v^fw^m i l r(n + 1/2)(2/ ° (n) - '■*<»» (91) 

for closed shell nuclei. The coefficients Io(n) and I a p(n) [24] are given in ap- 
pendix B. 

5. 3 Weight and Pseudo- Weight Functions for Closed Shell Nuclei 

The various ingredients occurring in the calculation of the weight function 
W^(z,cu) have been defined in terms of the spin-isospin and the eigenstates 
occupied in the state determinant D^^{x) describing the ground harmonic 
defining the state of the nucleus. The (a\n,£\fl) coefficients are known ana- 
lytically. Therefore, it is only a matter of computer programming to get the 
the needed W^(z,u) functions. Nevertheless, analytical calculations can be 
pushed further when we have to deal with closed shell or sub-shell nuclei by 
using the summation over the spherical harmonics. The analytical expressions 
are interesting because the nucleus is constituted by a core of nucleons in 
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closed shells and by other nucleons outside the core. The core is spherical 
with a total angular momentum J = 0. The particular state of the nucleus is 
defined by the configuration of the last open shell. The analytical derivation 
of the coefficients for closed shells are given in appendix B. 

Since closed shell nuclei are spherical only terms with £ = appear in (61) 
and (62). A shell A a is defined by the value A a = 2n a + £ a in terms of the 
HO quantum numbers of the occupied state. The contribution to the direct 
term in Eqs. (83-86) for each spin-isospin state when all the shells for which 
A a = 2n a + £ a are filled from A a = to the last shell A a = £ m is [24] 

n\ ^ \pj V P + 3 / V n + 3 - Pj 



By taking the four spin-isospin states and the normalization (87) into account, 
the contribution of the direct term in Eqs. (83-86) for closed shell N = Z nuclei 
is [24] 



([L m ]\u,0y d = 2I (£ rn ,n) 4<5 0e + 2(5 ae + 5 Te ) + 5 {ar)e 



(93) 



where e = 0, a, r, or. 

The contribution coming for the exchange term, denoted by I a p, cannot be 
obtained analytically. However, the procedure for its computation is explained 
in the appendix. For closed shells with N = Z spherical nuclei one has the 
relation 



([L m ]\u, Of = 2 j/ (C, v) [45 e + 2(5 ae + 5 re ) + 5 {ar)t 

- Iaf3^m, V) S 0e + 2(5 ae + 5 Te ) + A5 (f7T)e 



(94) 



for the normalization (87). 



For the Coulomb case, where only the proton states are taken into account in 
the calculation for I and I a p, we have the relation 



([L]|z/,A) c = 2/ (z/)-/ a/3 (z/) 



(95) 



A more compact expression for the weight functions W^ e \z) can be obtained 
by using the Rodrigues' formula (Ref. [23], Vol II p. 169) 



{-2) n n\P^\x) = (1 - x)~ a {l + x)- b D n x (1 - x) a+n {\ + x) 



\b+n 



(96) 
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where denotes the operator D™ = d n / dx n , a = L m + (D — 5)/2 — n, and 
b = 1/2. Then, for £ = Eq. (86) becomes, after integrating Eq. (62) over u, 



(e) _ 2-(^+ s / 2 - 2 ) T(L m + D/2) 

w iL m ]y z ) - - 



v^F T(L m + (D — 3)/2 



21 - n 



n=0 



Alternatively we may write 



( j v^V2j r(L m + (£> — 3)/2) 2j [L + Z) 



2£„ 



X 



]T(-2) n n! ([L m ]|n,0) (£) (l - z f^-n p L m +(D-^)/2-n^^ (g8) 



n=0 



The weight function can thus be written as a product 



(99) 



where now a = L m + (D — 5)/2 — 2£ m , (3=1/2, £ rn refers to the sum 2n + £ 
of radial and orbital quantum numbers in the last shell, and p^ m ^(z) is a 
polynomial of degree 2£ m , 



pfj{z) = C« £(-2)"n! {{L m ]\n^\l - z f^-n p L mH D-,)/2-n, P{z) (10Q) 



2£r, 



n=0 



The normalization of (97) is given by the first term, i.e for n — 0, 



+i 

/ W^)dz = ([L m ]|0,0) 

-i 



(0 



(101) 



where we used the fact that for n > 
+i 



J D n z [(l-^) a (l + ^) n+1 / 2 
-l 



dz = 0. 



The p^]!- 2 ); being a polynomial of degree 2£ m , can be expressed in terms of 
its roots 



(102) 



71=1 
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which can simplify the numerical calculations. 

A similar expression can be obtained in terms of the variable X = (l + z)/2 = 
rfj/r 2 , = (1,2), for < X < 1. By substitution one gets 



T(L m + D/2) 



v^r(L m + (L>-3)/2) 
x £([L m ]|ra, 0)< e > D n x {{\ - x) L ™ + ( D -VI 2 X n+l ' 2 ) dX (103) 



Analytical formulas for the hyperradial part of the two-body potentials written 
in terms of X can be obtained from (103) since 



J f(X)D x [(1 - X) a X n+l/2 ] dX = (-1)" J {I- x) a X n+l/2 dX (104) 



The integral can be obtained analytically when f(X) is either a polynomial 
or an exponential or the product of both. 



6 The Projection Function 



We have seen that in order to have an equation for the calculation of the 
two-body correlations, using Eq. (24) one must extract from the amplitude 
F(rke, r) where the pair particles (k, £) is not the reference pair (i, j), the part 
of F(rk£,r) which, nevertheless, depends on r^. There are two cases: Either 
one of the (k,£) is either % or j, in which case we have a connected pairs like 
(i, k) or (j, k) where k ^ i or j; or k and £ are neither % nor j and we have 
disconnected pairs. The procedure is the following i) Expand F(r^, r) in terms 
of HPs which are polynomials in r\ t and r 2 (which contains r 2 j). ii) Write r\ t 
in terms of the kinematical rotation vector (A. 3) [6,25]. iii) Extract from each 
polynomial the part which depends on r 2 - and ignore the residual part which 
contains many body correlations [8]. 

Practically since the two-body amplitude in Eq. (34) is expressed in terms 
of the z = cos 20 = 2r^/r 2 — 1 and r 2 (for = (1,2)), the expansion is 
done for any pair in terms of Potential Harmonics (PH) for pairs in S-state. 
They have been designed in order to provide a complete expansion basis for 
pairs of particles in S'-state. For bosons in ground state the PH are the Jacobi 
polynomials b ^ 2,1 ^ 2 (z) associated with the weight function Wq{z) (see Eq. 
(A.8)). 

Let z = cos 20 = 2r 2 j/r 2 — 1 for the reference pair (i — 1, and j = 2 in Eq. 
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(27)) and zu = 2r|^/r 2 — 1 for another pair. The two-body amplitude can be 
expanded as 

P(zkt,r)= E U/(l- Z 0°(l + ^WW W 1 r)dz'(lO5) 



where P^ ,l3 (z) are Jacobi polynomials associated with the weight function 
(l-z) a {l + zf with a = (D-5)/2 and /3 = 1/2. The /i^ is the normalization 
constant for the Jacobi polynomials 

h *J> 2^Y{K + a + l)T{K + (3+l) 
K ~V(2K + a + p + l)K\T(K + a + (3 + 1) 1 1 



The expression in the braces in (105) is the 5-function S(z' — z M )- 

It was shown in Ref. [6] that for the connected pairs where z c = z^ or z c = Zjk, 
k 7^ i or j, and for equal mass particles P^ ,l3 (z c ) can be separated into two 
terms 



P k\zc) = * i, 'J PS^z) + other terms (107) 



and for disconnected pairs Zd = Zke, k,£ 7^ 



k P (^) = P le\ II pk^z) + other terms (108) 



Since P£ (z) is a HH we have Ar 2K P^ (z) = whatever the pair is. There- 
fore, A[r 2K x {other terms}] = and the residual part is also a HH orthogonal 
to P%- (z). Since only Pr (z) contains two-body correlations, the other terms 
are related to many-body correlations that we neglect. Finally, the projection 
function for pairs in ^-states is the sum of the projection functions of the 
2(A — 2) connected pairs and the (A — 2) (A — 3)/2 disconnected pairs. Thus 
we may write 

+1 

V°P(z, r) = J /[o] (z, z')P(z', r) dz' (109) 

-1 



with 
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oo npa,/3/ i /n\ , A-3 pa,P ( i\ 

f [0] (z,z) = (l-z) + -j^-t 

xP$ /, (z)P2 P V)/h& (110) 
An analytical expression of the projection function has been derived [8,13]. 

The projection of a potential harmonic for a pair (k, £) on a reference pair 
both in 5-state, is given by [6] 

^25) a = {D _ mfJ = 1/2 
P a /(1) ' 1 

where 5 = 27r/3 for connected pairs and 5 = 7r/2 for disconnected pairs. The 
projection function in terms of the angular parameter 5 is for one pair 



uz,z',5) = (i-mi+zr t ? K^Zl 5) * a KW a KV)ih a / (in) 

K=Q 



The sum over the series, Eq. (Ill), can be carried out analytically [8,13] 
providing the projection function in terms of cos 5 which is related to the 
choice of the pair to be projected where 5 = 2n/3 or 5 = n/2 for connected 
or disconnected pairs of equal mass particle. In terms of cos0 = r 12 /r = 



\J(l + z)/2 for the reference pair 



+i 

J f [Q] (z,z',cos25)P(z',r) dz' 



4 r(A + l/2) 1 
y/n T(A) sin 25 sin 20 



x 



1 


2(A-1) b 

J 


sin 5 sin (f> 


a 



(112) 



where a = cos(0 + 5), b = cos(0 — 5), and A = a + 1/2 = D/2 — 2. For 
disconnected pairs since sin 25 — > we have to take the limit of Eq. (112) for 
5 — > 7r/2 the result being 



/ f [0] (z, z>, -l)P(z>, r) dz' = ^(^('V " z) l,2 - X 

—z 

x | [-(z + /)] A - 2 (lW) 1/2 P(*>)d*' (113) 
-l 

For fermions and in particular for nuclei one follows the same procedure as for 
bosons in S-state. One starts from the weight function W[ Lm \(z) to find the 
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associated polynomials P^ m ' (z) which fulfill the normalization condition 

+1 

/ P% m] (z)P%r ] (z)W [Lm] (z) dz = 5 KK , (114) 



For bosons in ground states we have seen that they are the normalized Ja- 

cobi polynomials P K ^ \z) / ' \]~h^f '. From P^ m \z) one generates the projection 
function 



f [Lm] (z,z') = (A-2)Y: 



2P^ ml (-|) + ^Px m] (-l) 



xP^(z)P^(z>)W [Lm] (z>) 



K=0 



p^] (l) 



(115) 



where, whatever the normalization of W[ Lm ](z) is, the polynomials P^ m \z) are 
normalized according to (114). 

Let be Z\ • • • , z 2 e m the zeros of P[L m ]( z )- According to the Christoffel's formula 
, the polynomial Pk^z) is given by the determinant 



PT\z) = N [ £ ml I P[L m] {z) 



> PK+n( Z ) 

a,0 



Pr ( z n), PK+l( z n), 



P 
i 1 j 



K+n\ Z n) 



(116) 



where iVj^ m ' is a normalization constant fixed by Eq. (114) and P^ ,/3 (z) are 
Jacobi polynomials. Obviously for bosons in ground state where L m = one 
recovers the Jacobi polynomials for a — (D — 5)/2 quoted in Eq. (105). 

The expression of P^ m \z) in terms of a determinant, Eq. (116), presupposes 
that the roots of the polynomial P[L m ]( z ) ( e = 0) are known. An alternative 
way to obtain P^ m ' (z) is to use the moments formula [23] 



j [in 

K 



( z ) 



1, (l + z)/2, 
Co, C\, 

Ck-i Ck, 



,((l + *)/2) 
, Ck 
j Ck+i 



K 



c 



2K-1 



(117) 
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where the elements Ck are given by 

'1 + z 



+ 1 ■ - K 



C K 

-i 



/(-^) ( U8 ) 



and Cj^ m ' is the normalization constant. 



7 The Effective Nuclear Potential 



The Eq. (34) for IDEA has been obtained from (27) by multiplying at left by 
D* Lm ^ (Q) and integrating over all angular coordinates Qn-i associated with the 

Jacobi coordinates £j for i < N, (£jv — ^12) and then by dividing by the weight 
function W[L m j(z). Eq. (34) thus obtained is valid for Wigner-type potentials, 
i.e, without taking into account the exchange operators. In general, however, 
nuclear potentials include exchange operators that give rise to pseudo-weight 
functions as well. Thus, the right hand side of Eq. (27), after integration 
over dfijv_i gives E e Ml u)W(£ n ) with £ N = r 12 = rJ(l + z)/2. For 



spherical nuclei the weight function does not contain the angular coordinates 
u of £n and we may define 



W \U*)= w ilUM^ ( 119 ) 



and the effective potential in the right hand side of Eq. (34) reads 

v(t N ) = v«(t N )= £ (120) 

e=0,a,T,(7T 



while the effective hypercentral potential in Eq. (34) is 

+1 +1 
V [Lm] (r) = J W^ m] (z)V cS (U) dz = J £ W$Jz)VM(£ N ) dz (121) 
-1 -1 ( e ) 



The nuclear potentials are given in terms of the triplet or singlet for even or odd 
states. The following relations hold between the V^ e \ where e = 0, a, r, or, 
and the V 3+ , V 1+ , V 3 ', and the V 1 - potentials 



V 1+ + V 3+ + V 1 - + V 3 -] (122) 
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1 r 


v a = 






4 




1 r 


V T = 




4 




1 r 


V UT = 





-V 1+ + v 3+ - v l ~ + v 3 ~ 

yl+ _ y3+ _ yl- + y3 

-y l + _ \/ 3+ + V 1 - + V 3 



(123) 
(124) 
(125) 



In terms of Wigner, Bartlett, Heisenberg and Majorana potentials, the follow- 
ing relations hold 



yW = yO^ yB = yc ^ yH = _yr ^ yM = _yc 



(126) 



For the N = Z nuclei constructed from spin and isospin saturated HO states, 
like the a-particle, the relation (94) still is valid and the effective potentials 
associated with the direct and exchange terms I and I a p are 



v Io = \ {v 1+ + v*+) + \v^- + 9 -v 3 - 

for Iq and 

v Ig = -- (v 1+ + v 3+ ) + V- + -V 3 



(127) 



(128) 



for I a p. 

The effective potential can be alternatively defined in terms of the polynomials 



(129) 



with a = L m + (D — 5)/2 and 



2(>„ 



Pitt*) = E^Ia,(n)(z - l)^-npa-n,l/2 {z) 
n 

the result being 



(E) 



V e s{rij) 



Ml1]( z ) ~ P[Z] 



(E) 



(130) 



(131) 



where Vd = 2V/ and Ve = 2V/ a/3 and where the polynomial occurring in the 
weight function is 



(132) 



38 



It is interesting to note that the one pion exchange potential (OPEP) defined 
by 

Vopep (Uj ) = (a i -a j ){T i -T j )Y (r i:j ) (133) 
where 

y(rij) = v — 

for 3Vo ~ lOMeV and /i ~ 0.7fm _1 does not contribute to the direct term. 
Indeed, since 



^opep — K^pep — 3Y(rij) 
V& EP = 9Y(r ij ) 

^OPEP = Y( r ij) 

the contribution in the direct term disappears and amounts to lSYfcj) in the 
exchange term, Eq. (128). 

Since the contribution to the effective potential of the exchange term decreases 
rapidly for increasing A, J (0) = (A/A) 2 , while I aP (0) = A/A the contribution 
of the OPEP, which is the dominant term in the description of the long range 
part of the nucleon-nucleon potential, fades away for large nuclei. 



7.1 Spurious component 



In the Schrodinger equation, Eq. (25), the sum over all pairs of the two-body 
potential contributes while only one component V(r,ij) associated with the 
reference pair appears in the amplitude of Eq. (24). The potential can be 
expanded on the complete potential basis associated with the weight function. 
The polynomials of degree one 

where iV is the normalization constant given by 

^/(^-j^V^. 

-1 

is independent of [L m ] and disappears when the sum is taken over all pairs, 
J2i,j>i r ij/ r2 — A/2. Therefore one term occurring in the expansion of the 
potential disappears in the Schrodinger equation. 
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The occurrence of this spurious component in the potential has been discussed 
in Ref. [18,27]. This component contributes in the SIDE, i.e when Vjx, m ](r) in 
Eq. (34) is set to zero since in the right hand side P(z,r) is an amplitude 
where the reference pair is in 5-state while V(rij)V° operates as an S- 
state projected potential that vanishes when in the other amplitudes expanded 
in terms of the Jacobi coordinates, Eq. (A.l), the reference pair is not in S- 
state [26]. 

For three-bodies this equation is known as the Faddeev equation for ^-states 
projected potentials. In the IDEA we assume that the potential is local and 
we isolate the hypercentral part of the potential, Vjx, m ](r), leaving only the 
residual potential to operate, in the r.h.s of Eq. (34), on pairs in ^-states. 
When the sum in Eq. (24) is performed over all pairs (i, j), one obtains the 
Schrodinger equation Eq. (25). 

Since the component P\{z) in the expansion of V{rij) disappears in the sum 
J2ij>iV( r ij) occurring in the Schrodinger equation, it must be canceled in 
V(rij) for local potentials. Therefore, in order to avoid taking spurious part 
of the potential operating in IDEA, one should consider instead 



V(ry)-K(r)Pi(z) 



(134) 



where 



+i 



Vi(r) 



J V(r^(l + z)/2)P 1 (z)W [Lrn] (z)dz 



(135) 



-i 



Since 




and for the reference pair 



V°P 1 (z) = P 1 (z), 



the projection on this pair is 



V° ]T Pi(z M ) = -Pi(z) 



e,k>e 



which is in agreement with the coefficient 



(A -2) 2Px(-l/2) + 



A-3 



Px(-l) /Pi(l) = -1 



2 



occurring in the projection function, Eq. (110). 
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8 Spin-isospin Exchange Generated Elements of the State 



Once the problem for the closed shell nuclei in ground state has been solved, 
with the inclusion of two-body correlations, one may argue that at the level of 
accuracy where many-body correlations are neglected, the solution obtained is 
complete. However, the exchange operators occurring in the potential generate 
new states of grand orbital L m + 2 which can not be reached by solving the 
single IDEA equation, Eq. (34). 

To see how to include these states, let us begin with the a-particle where 
the problem is well known (the same situation holds for three nucleons in 
ground state). For three- and four-body systems, besides the space symmetric 
state associated with the spin-isospin antisymmetric state, the spin and isospin 
exchanged operators can generate also mixed symmetry states which must be 
taken into consideration in constructing the fully antisymmetric state. In such 
a case we have a sum of space and spin-isospin states products each state 
of the product having a definite symmetry in the exchange of two selected 
particles. To understand how such states are generated let 

D a = \\ap j3p an (3n\\ (136) 

be the spin and isospin antisymmetric Slater determinant associated with the 
fully-symmetric space state of 4 He in ground state, where ct(fi) denotes the 
spin- up (down) while p denotes the proton and n the neutron. The ground 
polynomial, from which the expansion of the 4 He in ground state starts, is 
DoXm where Yp, is a constant, in fact, a HH of order zero, i.e, with L m = 0. 
Applying, for example, the Bartlett spin-exchange operator on D a one 
obtains either 1 when the two spins are the same or -1 when the exchanged 
spins between two identical particles (protons or neutrons) are opposite, or 
when the exchange is between two different spins belonging to two different 
particles. Indeed, in the last case the determinant has two pairs of identical 
columns and thus disappears, i.e, one gets that J2ij>iPijD a = 0. Similarly, if 
the isospin is exchange between either the first and the last column or between 
the second and third column for which C<* + C/3 — 0, one gets a determinant 
with two pairs of identical columns and thus also disappears. 

If now one considers the scalar operator J2ij>i r ijPij t ne exchange between the 
same columns generates new HPs, namely, 

Pipi aiini Q!iniXi|| (137) 



and 




(138) 
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where in the expansion of the determinant with respect to the first two rows 
the scalar product X\ ■ x 2 appears. These two determinants are HP of degree 
two which cannot be generated by the Wigner or Majorana potentials and 
should be included in the description of 4 He. These are the so-called mixed 
symmetry states which are coupled to the fully symmetric S'-state through the 
spin-exchange (Bartlett) potential. 

Apart from the scalar spin-exchange operator rf^-P? one can apply also the 
scalar isospin-exchange operator rfjP[j. The generated determinants 

\\anx j3p an j3px\\ (139) 



and 



\ap f3nx apx (3n\ 



(140) 



are HPs of opposite sign to those generated by the spin-exchange operator, i. e 
r*P£D a = -r*I*D a , (141) 



and thus one can consider only HPs generated by the spin-exchange operator 
and coupled to the space symmetric state D a through the potential 

(V^P* + VT{ rij )Pl 3 ) D a = \ (y 3+ (r tj ) - V 1+ (r, tl )) P°D a (142) 

where V^ 3+ ^ and V^ 1+ ^ are the triplet even and singlet even potentials respec- 
tively. It is clear that the HPs generated by the operator rfjP^D a cannot be 
neglected. They constitute the mixed symmetry component of the 4 He wave 
function coupled, through the spin- and isospin-exchange operators, to the 
space symmetric component of the wave function. 

Similar components are generated by the exchange operators in nuclei. They 
should be taken into account explicitly as they bring a non-negligible contri- 
bution to the binding energy. The potentials associated with the spin- or the 
isospin-exchange operators depend on the difference between the triplet V 3+ 
and singlet V l+ even potentials which, for realistic potentials, is rather weak. 
Indeed, the contribution to the binding energy of the mixed symmetry state 
in three- and four-body bound states does not exceed a few MeV. In contrast 
to the Bartlett potentials, however, the Majorana potential depends on the 
sum of these potentials which is large. It is responsible, together with the 
Wigner potential, to nearly all the binding energy in few-body systems. The 
spin-isospin exchange operator does not generate mixed symmetry states in 
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the tri-nucleon or the alpha-particle system but in nuclei it does and therefore 
it can not be neglected. 

Let us note by D a \n, £, m > the part of a system of particles where four 
nucleons are in the same space HO state \n,£,m > and the spin and isospin 
states are saturated as in D a (see Eq. (136)). For instance, for 16 the HO 
wave function would be written as 

Df L ° m] = D a \0, 0, 0)£> a |0, 1, -1) D a \0, 1, 0) D a \0, 1, +1) (143) 

where D a \n, £, m > means that each of the four spin-isospin states in D a is in 
the same space \n,£,m > state. Let us further consider the operator 

^ = E(l-^)^7 (144) 

i,3 



which exchanges the spin-isospin states Q and Q ^ Q. When the exchange is 
between spin-isospin states in different states \n, £, m > and \n', £', mf > it gen- 
erates a determinant with two identical pairs of columns and then disappears. 
Proceeding, as for for spin and isospin, the scalar operator 

E^-^a-w^; ( 145 ) 



where 

c f 1 when n a = n' a , £ a = id, m a = m' a 
1 otherwise 

applied to different space states \a > and \a' > generates HPs of degree L m + 2 
which should be included in the description of the wave function. They should 
much contribute to the binding energy of the nucleus since they are coupled 
to the main component of the wave function through the Majorana potential 
which contains the sum of the triplet and singlet even potentials. 

It has already been shown, by solving the 6 Li problem in the hypercentral 
L m -approximation [28], that the contribution to the binding energy of the 
mixed symmetry state generated by the Majorana exchange operator is of the 
order of a few MeV. For a HO Serber force which is fitted in order to give 
the experimental size and binding energy of 4 He, the eigen-energy calculated 
for 6 Li without the inclusion of the mixed symmetry state is above the one of 
4 He [28]. An exact calculation could not provide, for potentials vanishing at 
infinity, an eigen-energy above the one of 4 He since one should at least obtain 
the 4 He binding energy for 6 Li when the two nucleons in the p-shell are in zero 
energy scattering state. 

When the operator (144) operates on closed shell it generates two pairs of 
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identical columns. Therefore, the operator (145) in order to remove the can- 
cellation of the determinant must change the space state of the incriminated 
columns. It means that only the scalar product Xi ■ Xj in rf- must be taken 
into account. But the scalar product written as 

4tt +1 

Xi ■ xj = XiXjP^cos <p) = x iXj — ]T Y^iu^iuj) , (146) 

m=-l 

where tp is the angle between the directions uji and ujj and Pi (cosy?) is a 
Legendre polynomial, is a harmonic polynomial of degree two as the product 
of two harmonic polynomials of degree one in the variable Xi and Xj. 

When the operator 

qr = £ % ■ % [i - %cj(i - <W^7 (147) 

i,j>i 

is applied to the state -D^j (#) it generates a sum of new HO states (£) 
which are HP of degree L m + 2. The new polynomials are indeed constructed 
according to the procedure described in Sect. II for HPs. Let us call 

Dl[ m+2] {x) = Df L %x) 

the sum of the HO Slater determinant generated by the Q UT operator. Each 
one is constituted by the original D^^x) Slater determinant where two HO 
individual states in the last occupied shell have been raised to the next shell. 
The -DjjP +2 ] (x) must be normalized and the weight and pseudo- weight func- 
tions and projection function must be calculated. Then the coupling between 
the Dyj^{x) and D^^ +2 ^(x) states generated by the spin-isospin exchange 
operator can be introduced through the exchange pseudo-weight function 

J D^ m+2] (x)P^Df L Z } (S)d^ N _^ 



9 Applications 

To solve the basic equation (34) for a specific nucleus and a specific filling of the 
HO shells, one has first to evaluate the coefficients ([L m ]\n, from which 
the weight functions (99) and the density functions Az), Eq. (100), or 
equivalently their roots (102) can be obtained. These coefficients are evaluated 
from Eqs. (83), (84), and (85) as well as from Eq. (89) when Coulomb forces 
are included. The results for the closed shell nuclei A = 16 and A = 40 and for 
the closed shell neutron systems N = 8 and iV = 20 are given in table 1 while 
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the corresponding roots for pf^Az) are given in table 2 for the 16 and 40 Ca 
nuclei and in table 3 for the closed shell 8-neutron and 20-neutron systems. 
It is noted here that only the relative normalization constants C^/C^ are 
required to construct the effective potential defined by 

e=0,a,T,<jT 

- E V^\r l3 ) P $Jz)/p$Jz) (148) 

e=0,a,r,<TT 

where the nuclear potentials are given by Eqs. (122-125) in terms of the 
singlet or triplet, even or odd potentials. The extracted ratios C^/C^ are 
also given in tables 2 and 3. It is further noted that the weight function for 
the 4 He system, is given by (99) with p^ e \z) = 1 for L rn = 0. 

The next task is the choice of the potential. Since we are testing here the 
suitability of our method in nuclear structure calculations, we employed three 
interactions with different characteristics and for which results by other com- 
peting methods are available. The first potential employed is the widely used 
in nuclear structure calculations Brink and Boeker Bl effective soft core po- 
tential [29] which is of a rather long range and of soft core. Although this 
potential is not realistic, since it does not fit the scattering N-N phase shifts, 
there exist in the literature a lot of results obtained with it and therefore its 
use is warranted for comparison purposes. For convenience, we recall here this 
potential which has only Wigner and Majorana components 

V w = 595.55 exp(-2.041r 2 ) - 72.212 exp(-0.512r 2 ) , 
V M = -206.04 exp(-2.041r 2 ) - 68.388 exp(-0.512r 2 ) . 

The second potential used is the Afnan and Tang [30] S3 potential. The original 
potential was adjusted to the static properties of 4 He nucleus. However, in 
order to extend its applicability to heavier nuclei, Guardiola and collaborators 
[31] added a repulsive part in the singlet- and triplet-odd components, the 
modified potential thus obtained (known as MS3 potential) being 

V w = -5.75 exp(-0.4r 2 ) - 10.75 exp(-0.6r 2 ) 

- 41.5 exp(-0.8r 2 ) - 81.675 exp(-1.05r 2 ) + 1000 exp(-3r 2 ) 

V M = 5.75 exp(-0.4r 2 ) + 10.75 exp(-0.6r 2 ) + 41.5 exp(-0.8r 2 ) 

V B = -V H = 5.75exp(-0.4r 2 ) - 10.75 exp(-0.6r 2 ) 
+ 41.5exp(-0.8r 2 ) - 81.675 exp(-1.05r 2 ) 

This potential has a very strong repulsive core which give rise to strong two- 
body correlations. As a third potential we employed the more realistic soft 
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core interaction of Gogny, Pires, and de Tourreil (GPDT) [32] which was also 
used in the past in nuclear structure calculations. 

Using the aforementioned potentials we firstly investigated whether 8 and 20 
neutrons, forming a closed shell system, can sustain a bound state. It was 
found that, at least with the potentials employed, no bound state exists. Even 
with the soft Bl potential a bound state can only be generated by modifying 
the potential by more than 20%. 

The ground state energies E g and root mean square radii r rms for the closed 
shell nuclei 4 He, 16 0, and 40 Ca are given in tables 4, 5, and 6. The results were 
obtained in the extreme and uncoupled adiabatic approximations and are com- 
pared with other results obtained with cluster expansion method (FAHT) [31], 
Bruekner-Hartree-Fock type (BHF) [31], Fermi-Hyper- Netted-Chains (FNHC) 
[33], variational Monte Carlo (VMC) [34,35] as well with results obtained with 
hyperspherical harmonics expansion methods (HHE) [10]. We refer also to the 
work of Guardiola et al. [36] where more relevant results are compiled using 
various techniques. The importance of the correlations, stemming mainly from 
the short range repulsion of the two-body force, is also inferred from these ta- 
bles by comparing the IDEA results with those of hypercentral approximation. 

It should be noted that the results obtained by our method do not involve any 
adjustable parameter unlike, for example, the FHNC method whose results 
depend on the Jastrow ansatz, and in general on the model wave function 
employed. This can be seen, for instance, in table 5 where the use of different 
ansatzs in the FHNC resulted in ~13 MeV difference in the ground state energy 
with Bl potential. 

The inclusion of the Coulomb potential is of utmost importance when one 
is dealing with nuclei involved in reactions of astrophysical interest. Its in- 
clusion, however, with the exact integral, differential, or integro-differential 
methods employed in Few-Body calculations is non-trivial and thus it is usu- 
ally omitted. However, the incorporation of Coulomb forces in our formalism 
is straightforward. The extra repulsion generated by the Coulomb potential is 
also given in tables 4, 5, and 6. As expected, the smaller the rms radius, the 
higher the eigen-energy is. 

We calculated the EAA and UAA giving respectively a lower- and upper- 
bound to the eigen energy. The difference decreases from ~ 0.9 MeV to ~ 
0.1 MeV for A growing from A=4 to A=40. By taking the average for the 
eigen energy, the difference with respect to the exact value is less than 0.5 MeV 
for 4 He and becomes negligible for increasing A. The values obtained are in 
agreement with the spectrum of those obtained by other methods except for 
40 Ca with the MS3 potential where we got a significant lower binding but 
unfortunately variational values are not available for comparison. We notice 
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that the strongest the core of the potential, the larger is the increase of bind- 
ing energy brought by the correlations. Since most realistic potentials have 
strong repulsive core, we expect that the binding in nuclei originates from the 
correlations and that the effect of the hypercentral potential is to balance the 
kinetic energy only. 

As an example for the applicability of our formalism to open shell nuclei we 
consider the 10 B nucleus. The ground state of this nucleus is known to be 
J w = 3 + , T = 0. Two likely ground configurations with this state are the 

ipi = (4si/ 2 )(4pi/ 2 )(2p3 /2 ) 

configuration where the p±/ 2 nucleons have the quantum numbers J = and 
T = and the 

^2 = (4si/ 2 )(6p 3 / 2 ) 

configuration where the two holes in the, otherwise, full subshell are also 

coupled to J w = 3 + , T = 0. The form of the density function P[L m ]( z ) i n this 
case is given by [24,5,37], 

= E (Wl^) (f) (l -z) L —(l +z) n (149) 

n=0 



where the expansion coefficients are given by 

= r ( I. + J / r 2 ( ^ + -3/Sr(n-3/2) < ■ (150) 



The A € n are calculated from the shell structure of the nuclear state under 
discussion by means of the Talmi-Moshinsky or Gogny coefficients [38,39]. 
The explicit values for these coefficients are given in table 7 while the roots 
and the relative normalization constants for the polynomial p^ m ] are given in 
table 8. The corresponding binding energy results for the two configurations 
are given in table 9. 



10 Conclusions 



When the structure of nuclei has been identified as associated to the HO 
quantum numbers corrected by the effect of a strong spin-orbit force, it was 
taken for granted that nuclei can be described as an IPM. In the present 
work we have shown that the HO is both an IPM and a collective model 
and that the HO potential is a hypercentral potential invariant by rotation in 
the D-dimensional space spanned by the particle coordinates in the center of 
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mass frame. This leads to another interpretation where the most significant 
part of the interaction is hypercentral in such a way that in the ground state 
the repulsive centrifugal barrier must be minimum. This interpretation also 
leads to the appearance of the so-called magic numbers in nuclei whatever the 
nuclear potential is. 

We have further shown that the next improvement after the HCA, where 
the wave function is described by the product of an antisymmetric harmonic 
polynomial and a function of the hyperradius, is obtained by introducing the 
two-body correlations generated by the nuclear two-body potential. This re- 
sults to an integro-differential equation (IDEA) for the A-body bound system. 
In order to test the quality of this approximation where only one harmonic 
polynomial of minimal degree and two-body correlations are taken into ac- 
count in the wave function, we compared the binding energies obtained with 
the IDEA to those computed for central forces by other methods. We found 
values in agreement with the spectrum of those available in the literature, the 
largest differences being in the case of 40 Ca with S3 potential which has a 
strong repulsive core. Unfortunately, no variational results are available in the 
latter case. 

Our formalism is based of course on the assumption that the NN force is the 
dominant one inside nuclei. It is well known, however, that three-body forces 
are also essential in the description of nuclear dynamics and their introduc- 
tion might explain at least a part of the underbinding of nucleons in nuclei 
interacting by two-body forces only. Their employment in our basic equation 
(34) is straightforward when they are symmetrical and separable, 

W(x)= £ V^(r t3 )V^(r jk ), 

i<j<k<A 

as described in Ref. [13,16]. 

Up until now it was possible to solve three- and four-body problems with 
accuracy by using the Faddeev and Faddeev-Yakubovsky equations respec- 
tively. For more nucleons it is customary to assume that nucleons in nuclei 
form clusters, and then to construct effective inter-cluster interactions for the 
sake of solving the problem as a three- or four-body. The construction of the 
effective interactions, however, is a formidable task with all sort of ambiguities 
naturally creeping in due mainly to insufficient scattering data required in the 
construction of the nuclear forces involved. 

The existence or not of bound and of resonance states in three- and four- 
neutron has been the subject of numerous theoretical and experimental works 
[40,41]. Going beyond the four-neutron system [42] is also important as nuclear 
matter exhibits a much richer phase structure than light nuclei because of 
stronger correlations between nucleons. However, our results for the 8- and 
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20-neutron systems indicate that realistic nucleon-nucleon central forces alone 
are not sufficient to generate closed shell bound systems. 



The huge differences between the results obtained with the three different 
two-body forces used, indicates that one has to employ a realistic two-nucleon 
potential supplemented with a three-body force that accounts for the missing 
binding energy. The use of effective forces, either at the two-body level, or 
cluster-cluster potentials, or optical potentials in processes where the details 
of the interior part of the wave function plays an important role, could provide 
misleading results. Such processes are, for example, the electromagnetic reac- 
tions, in which the transitions are guided by strict selection rules and depend 
on the details of the wave functions in each channel and their characteristics 
in the interior region. Therefore, the inclusion of the underlying correlations 
in an unambiguous way is of crucial importance. 



A Jacobi Coordinates 



Let us define the Jacobi coordinates for equal mass particles 



£at_i = ^(£3 - X 3 ) 



2i 

Cn-i+i = y ^ j [Xi+i - Xi) 



— ~(xi+i — Xi+i) 




(A.1) 

with N = A — 1, Xi = h J2p=i x p , and 



N 



r 2 = EC 2 = 2E(^ - X) 2 = \ E % - , (A.2) 

i=l i=l i,j>i 



and where X = X A corresponds to the center of mass. The normalization is 
such that for % = N one gets the vector f 12 . 
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For describing the coordinates one introduces the kinematical rotation vector 
[6,25] 

— * — * — * 

f(5, if ) = cos 5£n + sin5cos(/9^7v_i + (terms in j < N — 1) (A. 3) 

— # — # , — 

where 5 and y? are angular parameters, ^at = — Xj, and £at_i = \/3{xk — 
l/3(xi + xj + Xk)). For 5 = 0, r(0) = Xi — Xj is the reference pair, for 5 = 27r/3, 
99 = r(27r/3) = x^ — X; L is a connected pair, while for 5 = ir/2, r{ji/2) is a 
disconnected pair. 

One may introduce also the hyperspherical coordinates of Zernike and Brinkman 
[43] 

£ = 6v = f cos0Ar , p = rsin0 A r (A. 4) 



with 



£ 2 + p 2 = r 2 , P 2 =E^- (A.5) 

J=2 



For the other Jacobi coordinates we have 
£jv-i = rsin0Arcos0Ar_i 



£j = r sin • • • sin <j> j+1 cos 0^ , 

(A.6) 



£ 2 = r sin 0jv • • • sin 3 cos 2 , 
£i = r sin 0jv • • • sin 2 



where 0i = 0. It is clear that in the Zernike-Brinkman system of coordinates 
one has, apart from the hyperradius r, N — 1 angles 0j and 2N angles u>i, the 
volume element being 

dtt i+1 = dQi sin (/ ~ 4) i+1 cos 2 <f> i+1 d0 m du i+1 (A. 7) 



where I = 3(i + 1). The Q coordinates are separated into two parts, first the 
z = cos 20, (0 = 0jv), and uj = uj n , the angular coordinates of £jv, and second 
the fijv-i for the other hyperspherical coordinates (0j,cjj), i < N where Ui 
are the angular coordinates of Let df2Ar_i be the surface element of the 
unit hypersphere p — 1 in the D — 3 dimensional space spanned by the Jacobi 
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coordinates £j, % < N. The volume element in the D-dimensional space is then 
given by 



d^f = r"- i dr dQ 
dfi = (sin 4>) D ~ A cos 2 d0 duo dVl N -i 

= 2^(1 - ^ 5)/2 (l + ^) 1/2 d^ dc dlW 



= W n te) dz dw dfi 



where 



N-l 



COS (J) = — = - 
r r 



cos 20 = 2^f - 1 



(A.8) 



In general the surface element df2j is defined to be the part of the surface dfi 
which contains the coordinates ujj and Zj for j < i, i.e 



dfii = duji Y[ 2 Zj,2 Wj{zj)6.Zj^jjj 

3=2 



(A.9) 



where Wj(zj) is defined by 



it; 



i {z j ) = (l-z j )to-*V'(l + z j )W 



(A.10) 



Knowing the element dflj we may construct the element of dflj + i via 

dfij + i = dfij(sin0,,- +1 ) D_4 cos0 2 +1 d0j +1 dc<jj +1 (A. 11) 



where D = 3(j + 1). 

The kinetic energy operator T, for equal mass particles, is given by 

h 2 



n 2 A_1 



m 



3=1 



2M 



-V 2 



(A.12) 



We are interested to express the translationally invariant part of this equation 
in terms of the {u>i, Zi}. Due to the product structure of dVL it can be written 



as 



m 



d 2 3A-A0 L 2 (Q) 



dr 2 



dr 



(A.13) 
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where L 2 (Q) is the grand orbital operator. Let Lf(Qi) be the operator asso- 
ciated with the first vectors {£1, • ' ' > It can be written in terms of the one 
for the vectors {£i> " ' > 6-1} as follows 



L?(nO = £>(*) + + t^- ^(Hi-i) (A.14) 



where 



"W^^ 1 (A - 15) 



and 

L 2 (Q) = L^Q) . 
The / 2 (cUj) is the orbital operator normalized according to 

[F{u)+l(l+l)\Y lm (ui) = 0. (A.16) 



B Various Coefficients 



5.i T/ie ?>L- Coefficients 



The 3L-coefficients are defined by the expansion of a product of two HO 
eigenfunctions [44] 



7T 



n,£,m 



n a np n 

e a £/3 £ 



(B.l) 



where a and (3 stand for the HO quantum numbers n a , £ a , m a and np, £p, mp. 



By integrating over the x space for the HO quantum numbers a, (3 with 
n — 0, £ = 0, m = 0, and 



^,,0,0(5) = ^171 e-* 2 / 2 r °, 



(B.2) 



one finds the normalization 



n n 



np 0" 



= 5, 



aP 



(B.3) 
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Multiplying Eq. (B.l) by 

n!r(3/2) i 1/2 



_r(n + £ + 3/2) 
and integrating over x one gets [45] 



rr>)*X +1/ V) 



n a rip n 




_L £ P l 





n, 



1/2 



T(n a +£ a + 3/2) T(n p + ip + 3/2) T(n + £ + 3/2)_ 
x J x^^L^I\x 2 )L^I\x 2 )L^I\x 2 ) e~ x2 x 2 dx . (B.4) 



For n = £ = one finds the normalization (B.3). Using the expansion for the 
Laguerre polynomials for a and (5 gives for the integral in Eq. (B.4) 



1 = 



( _ 1)m w / Wq +£a + 1/2 \ / n/J + ^ + 1/2 \ 
m!m'! \ n a - m )\ np-m! ) 



This last integral is given, according to Eq. (58), by 

(-l) n /V 



Jx 2 ^L e n +1 / 2 (x 2 )x 2 dx 



2 \n 



\T(u + £ + 3/2) 



(B.5) 



where here v — m + w! + (£ a + ip — £)/2. 
The 3L-coefficients are thus given by 



n a np n 

£ a £/3 £ 



:-iy 



' 7T 



1/2 



X 



X 



n. 



np\ 



,1/2 



T(n a + £ a + 3/2) r(n^ + lp + 3/2) r(n + £ + 3/2) _ 

-l) mW Ai Q + 4 + l/2\ /% + ^ + l/2\ 
m\m'\ \ n n — m J 



E 



— m' 



x (™ + + {t l + C " - T(,n + m < + ^±1±1) (B.6) 
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We recall that (fy denotes the binomial coefficients, and that < m < n a , 
< m' < np, and n < m + m' + (£ a + £$-£) /2. 



B.2 The (a\n, £\(3) coefficients 



To obtain the coefficients (a\n,£\(3) which appear in the Fourier transform 
(74), we expand the plane wave 



oo +e 
e=o m=-e 



(B.7) 



where (k,u>k) and (r, ou) are the polar coordinates of k and f respectively and 



ji is the spherical Bessel function, je(p) = y 7r/2pJ^ + i/ 2 (p). Then, using the 
integral 



[(*) Lf+VV/fr r W/ A e^V 
J V o / v A;x 



6 3 y 2n - 



dx 



2 n+(^+3)/2 n | 



■e" y2/2 (B.8) 



where y = kb/y/2, one obtains 



{a\n,£\(3) = 



2tt 3 / 2 



nl/2 



22n+/ n |r( n + £ + 3/2) 



/ym a I ym* I ^ m /3 \ 
\ £ a \ £ I / 



n Q n/3 n 

4 ^ £ 



;b. 9 ) 



The (Y £ ™ a |Y £ rrt *|Y^™' 3 ) can be easily expressed in terms of the Clebsch-Gordan 
coefficients [46]. 

For n a = rip = 0, as for instance for the first p— and ci-shells, 



(0, £ a , m a \n, £, m\0, £p, mp) 



x 



2«+ £ / 2 r(n + £ + 3/2) 

\£ a + £/s- £)/2\ r((4 +£{, + £ + 3)/2) 



n ;[r(4 + 3/2)r(^ + 3/2)] 



{Y ^\Yr\Yr) (B.10) 



with m = 171/3 — iTia, and |£ Q — tp\ < 2n + £ < £ a + ^ The coefficient (B.9) 
is the product of two terms where only the first depends on the azimuthal 
numbers m a and mp but not on n a and np. 
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B.3 The Io(£ m ,n) coefficient 



For closed-shell nuclei the coefficient I (£ m ,n) can be calculated from (92). 
Table 10 gives the first 11 values of Iq for £ m — 0, • • • , 5. The rest can be easily 
calculated either from (92) or from the general relation [24] 



n\,ri2,ti 



^(a|ni,£i) 



(B.ll) 



with n — n\ + n 2 + £\- 



B.4 The I a p(£ mi n) coefficient 



They are given by 

l af >(t m ;n,0) = Yl £H^i|/3)(/3MiN 



«l,"2A a/3 

n—n 1 +n2+tl 



(B.12) 



where the sum over a and /3 means a sum over all H.O. states (each must 
be taken once only either for a or f3). The calculation of this coefficient is 
tedious. However, for closed-shell nuclei and for small values of n it reduces 
to the following analytical expressions: 



I a p(£m, 3) 

I a p{£mA) 
I a p{£m, 5) 
la/3 {£m, 6) 



3! 



Iap{£m, 2) 



1 (lm + 3)! = -A 

M 4 ' 



1 (C + 3)! _ A 
A\{£ m -1)\ 16" m ' 

1 (£rn + 3)\ 

5! (4,-2)!' 
1 5C-l(C + 3)! 
6! 6 (4 -2)!' 
1 £ m + l(£ m + 3)\ 
6! 12 (4,-3)!' 
94, 34. 32 



5!(C + 1) 
/ 1 \ 2 

4/?(4; 24) = yj, j" J j 



(B.13) 
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^(^mj 2£ m 1) 



1 



(£ m 2) ! £ m \ 



Iap(£m,n > 2£ m ) = 0. 
Values of these coefficients of up to £ m = 7 are given in Table 11 

B.5 Coefficients in the jj -coupling scheme 



For describing the state of the nucleus, the Slater determinant (2) has been 
filled with individual HO states each associated with its spin-isospin state. 
Instead, one can combine the spin s and angular momentum £ to generate the 

— * 

total angular momentum j of projection m and describe the individual states 
in terms of the isospin \t > and the \n,£,j,m > states. The total angular 
momentum is constructed as usual, 



(B.14) 



m,m s 



where Xm a is the spin variable with m s = ±1/2 and the bracket is a Clebsch- 
Gordan coefficient. 



The description of the individual states in terms of angular £ and total angular 
momentum j is well adapted to nuclei because the two-body spin-orbit opera- 
tor is actually a sum over individual spin-orbit operators (see Eq. (13). In the 
jj-coupling scheme we have to calculate the new coefficient (a\n,£\(3) where 
this time \a > stands for \n a , £ a , j a ,m a > and \/3 > for \np,£p,jp,mp >. The 
HO Slater determinant is now 



(B.15) 



where 



2n! 



b 3 T(n + £ + 3/2) 



1/2 



1^")lW V V) e-^ /2 (B.16) 



where y = x/b. 

The (a\n,£\/3) coefficients in jj-coupling must be constructed from the HO 
individual functions (B.16). One starts from the matrix elements where the 
integral is performed over doo 
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(L,ja,m a \Y e m *\e0,jp,mp) = ^=(-lH" 1/2 lijpijajf, 



1/2 



X 



-a ""[) 



J/3 3a 



/ 1 — mp m a m 



3(3 3a 1/2 



(B.17) 



in terms of the 3j and 6j coefficients. Here x — 2x + 1, The one-body co- 
efficient (a\n,£\(3) in jj-coupling is obtained by substitution of (B.17) for 
lJiy\Yr\YTf) ^ Eq. (B.9) 



(a\n,e\P) 



2tt 3 / 2 



1/2 



2 2n + i n\T{n + £ + 3/2)^ 
x (£ a ,j a ,m a \Yp*\ep,jfi,mfi) 



n a np n 

4 lp £ 



When i = then 

{£aJa,ma\Yo\£pJp,mp) 



Sail 



(B.18) 



(B.19) 



and 



(a\n,0\P) 



£) V4 [2 2 "n!r(n + 3/2) 1_1/2 



n Q n,3 n 
4 ip 



^a^^aj/3^m a m fl (B.20) 



In the direct term the sum over — j a < m a < j a for a complete subshell 
n a ,£ a ,j a gives 



£ (a|n, 0) = (2j a + 1) Q) V4 [2 2 "n!r(n + 3/2) 1 " 

m a =-ja 



;b.2i) 



In the jj-coupling scheme the Slater determinant is constructed according to 
(B.15). Instead of the four spin-isospin cases occurring in Eq. (70) we have 
to consider only two cases for dij(xi,x 2 ) in Eq. (63): Either the two nucleons 
are the same or we have to deal with the neutron-proton pair. In the first 
case P — P = 0. For the Wigner (e = 0) and isospin exchange e = r weight 
and pseudo-weight functions, the matrix element (B.18) must be used in Eqs. 
(79) and (81) to obtain the direct and exchange coefficients. The coefficients 
([L]|i/, A)( e ) are given by Eq. (83) for the weight function (e = 0) and by 
Eq. (85) for the isospin exchange (e = r) pseudo-weight function where t is 
substituted for (. 
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The spin exchange operator generates new HO states ^ a '{xi)ippi{xj) with 
PtjM^M^j) = E (a/,P\P? j \a,P)il> ,&)M*j) ( B - 22 ) 



3 a 1 ,3f} ' J 



where 



(a',f3'\P°\a,[3} = 



JaJpJa'J/3' 



(_l)J>+J> 



x (jaj/3m a mp | J M) (jadpim^mp, | J M) 



4 1/2 ja 
1/2 lp 3(3 

3a' 313 J 



(B.23) 



where j = 2j + 1, ot = \n a t£ a >j' a ,m a /}, f3' = \np'£p'j'p,mp'), the brackets are 
Clebsch-Gordan coefficients, while the brace is a 9j coefficient. 

When the states \a > and \(3 > are coupled to give a definite total angular 
momentum J of projection M, like for two nucleons outside a closed shell (e.g 
6 Li), the effect of the spin-isospin exchange operator is 



x ]T (-l)w+w 



1/2 



3a' dp' 



3 a j 13 



1/2 



1 m: 1/2 J a 

1/2 ^ j> 



M 



When the spin exchange occurs, the direct coefficient is given by Eq. (79) 
where 



(a,P\n 1 ,e 1 ,n2,e 2 \aP) {0) = 

J2 (<*', P'\P^a, P) (a|n 1; ^a^K £ 2 \P') (B.24) 

ja' ,3p' J 

where (a\n, £\ct') etc, given by (B.18), is substituted for (a|ni, £±)*(P\n2, £2)- 

The exchange coefficient is obtained by using ((3a\ni£\n2£2\ot(3)^ instead of 
(a.\n\£\\P)*(P\n2£2\ot) in Eq. (81). The spin (e = a) and isospin (e = ar) 
exchange coefficient ([L]|i/, A)^ and ([L] \v, A)^ are given respectively by 
(83) and (85) where t is substituted for (. 
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Table 1 

The coefficients (L m \n, l)^ for the closed shell nuclei A= 16 and A=40 and for the 
A=8 and A=20 closed shell neutron systems. 



e 


16Q 


40 Ca 


8n 


20n 


WIGNER 


120. 


780. 


28. 


190.0 




-64. 


-800. 


-16. 


-200.0 




6. 


270. 


1. 


65.0 






-38. 




-9.0 






1.5 




0.25 


BARTLETT 


48. 


360 


8. 


80.0 




-32. 


-400. 


-8 


-100.0 




0. 


120. 


-1. 


25.0 






-16. 




-3.0 











0.250 


HEISENBERG 


48. 


360. 


28. 


190.0 




-32. 


-400. 


-16. 


-200.0 




0. 


120. 


1. 


65.0 






-16. 




-9.0 











0.25 


MAJORANA 


0. 


120. 


8. 


80.0 




-16 . 


-200. 


-8 . 


-100.0 




-6. 


30. 


-1. 


25.0 






-2. 




-3.0 






-1.5 




0.25 


COULOMB 


28. 


190. 










-16. 


-200. 










1. 


65. 












-9. 










0.25 
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Table 2 



Roots and relative normalization constants C^/C^ for 16 and 40 Ca systems. 



C 


16Q 


40 Ca 




Roots 


c( £ )/c(°) 


Roots 


c«/c(°) 


WIGNER 


-1.2590158014 

-1.J.40ZOOZ40U 


1.0 


-1.1554853491 

-1.U10Z 1 UZ 1 4U 

(-0.9828722545, 0.0692920115) 
(-0.9828722545, -0.0692920115) 


1.0 


BARTLET 


-1.0 
1 n 

l.U 


-0.2085327360 


+1.0 
i n 

-l.U 

(-0.9723433872, 0.0576076251) 
(-0.9723433872, -0.0576076251) 


-0.1036460531 


HEISENBERG 


-1.0 
1 n 

l.U 


-0.2085327360 


+1.0 
i n 

-l.U 

( -0.9723433872,0.0576076251) 
( -0.9723433872,-0.0576076251) 


-0.1036460531 


MAJORANA 


-0.9283077625 
n (^971 dqi i n 


-1.5213318399 


-0.9879721459 

(-0.9604317388, 0.0448314378) 
(-0.9604317388, -0.0448314378) 


-1.2591151327 


COULOMB 


-1.9906725214 
-1.0483085670 


0.1319112107 


-1.4266427288 

-1.0070162824 
(-0.9781410398, 0.0634856828) 
(-0.9781410398, -0.0634856828) 


0.1493923245 
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Table 3 



Roots and relative normalization constants C^/C^ for 8n and 20n systems. 



C 


8n 


20n 




Roots 




Roots 




WIGNER 


-8.6346734966 

111 9^789^1 


1.0 


-2.1965713866 

1 01/1^0^9791 

(-0.9502197308,0.1301817431) 
(-0.9502197308,-0.1301817431) 


1.0 


BARTLET 


-0.9146656586 

n D997Q991 C\A 
-u.uzz / yzoiU4 


-10.3986013985 


-0.9870485512 

D 9/L9/L97D979 

(-0.9502197308, 0.1037682674) 
(-0.9502197308,-0.1037682674) 


-2.4763777976 


HEISENBERG 


-8.6346734966 
1 1 1 9^789^16 


1.0 


-2.1965713866 
1 D14^Qfi9791 

(-0.9502197308, 0.1301817431) 
(-0.9502197308, -0.1301817431) 


1.0 


MAJORANA 


-0.9146656586 
-0.0227923104 


-10.3986013985 


-0.9870485512 

-0.3424370373 
(-0.9270670083, 0.1037682674) 
(-0.9270670083, -0.1037682674) 


-2.4763777976 
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Table 4 

Ground state energies E g (in MeV) and root mean square radius r rms (in fm) for 
4 He. The IDEA results are those of the uncoupled adiabatic approximation (UAA). 
The results a in FHNC are with Gaussian and b with Euler correlations. Ec is the 
increase in the binding energy due to the Coulomb interaction. 



Potential 


HCA HCA+Coul EAA UAA E c 


Other Methods 


Bl[29] 


29.292 28.473 39.169 38.162 0.810 
(1 497 s ) (1 502") (1 420) (1 422) (1 425) 


FHNC/l a [33] 37.7 
FHNC ll b f33l 37 9 
FAHT/III [31] 36.6 
BHF [31] 36.9 
VMC [34,35] 36.4 


MS3[30,31] 


7.177 6.516 28.070 26.760 0.805 
(1.894) (1.910) (1.443) (1.454) (1.458) 


FHNC/l" [33] 24.7 
FAHT/III [31] 24.2 
BHF [31] 25.0 
VMC [34,35] 23.9-26.5 


GPDT[32] 


14.199 13.436 18.848 18.175 0.757 
(1.651) (1.661) (1.593) (1.603) (1.608) 
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Table 5 

Same as Table 4 but for 16 O 



Potential 


HCA HCA+Coul EAA UAA E c 


Other Methods 


Bl[29] 


106.529 93.135 164.777 164.332 13.727 
(2.604) (2.628) (2.559) (2.566) (2.587) 


FHNC/l a [33] 150.4 
FHNC/1 6 [33] 152.4 
FAHT/III [31] 163.7 
BHF [31] 163.7 
VMC [34,35] 150.9± 0.3 
IDEA[47] 165.2 
HHEflOl 152.10 

L J 

E c [10] 13.78 


MS3[30,31] 


12.745 1.443 103.261 102.793 13.843 
(3.079) (3.167) (2.539) (2.535) (2.563) 


FHNC/1 6 [33] 105.3 
FAHT/III [31] 107.7 
BHF [31] 118.6 
IDEA [47] 103.2 


GPDT[32] 


61.867 47.740 100.796 100.688 14.305 
(2.4791 (2.513) (2.484) (2.484) (2.528) 


HHE[10] 94.63 
E c [10] 14.71 
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Table 6 

Same as Table 4 but for 40 Ca 



Potential 


HCA 


HCA+Coul 


EAA 


UAA 


Ec 


Other Methods 


Bl[29] 

L J 


323.355 


250.844 


475.381 


475.345 


75.461 


FHNC/l a [33] 


471.0 




(3.341) 


(3.397) 


(3.251) 


(3.251) 


(3.300) 


FHNC/1 6 [33] 
FAHT/III [31] 
BHF [31] 


482.0 
478.0 
507.2 














VMC [34,35] 483.0± 0.4 














WFA [47] 


447.8 














IDEA [47] 


483.0 














HHE [10] 


468.14 














E c [10] 


76.0 


MS3[30,31] 

L > J 


48.419 


no bound 


259.447 


259.328 


72.149 


FHNC/1 6 [33] 


350.0 




(3.868) 




(3.375) 


(3.372) 


(3.457) 


FAHT/III [31] 
BHF [31] 
IDEA [47] 


335.6 
354.0 
272.54 


GPDT[32] 


255.952 


171.852 


376.331 


376.110 


86.31 


HHE [10] 


363.53 




(2.868) 


(2.951) 


(2.842) 


(2.842) 


(2.912) 


E c [10] 


86.0 
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Table 7 

The coefficients A e n , Eq. (150) for the two configurations considered for the 



nucleus. 



e 


(4a 1/2 )(4p 1/2 )(2p 3/2 ) 


(4s 1/2 )(6p 3/2 ) 




A £ 


A\ 


A e 




A\ 




WIGNER 


41/2 


19 


11/2 


41/2 


19 


11/2 


BARTLET 


1/2 


46/3 


1/2 


1/2 


40/3 


1/2 


HEISENBERG 


-1/2 


12 


-1/2 


-1/2 


12 


-1/2 


MAJORANA 


-171/8 


323/12 


-51/8 


-156/8 


284/12 


-36/8 
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Table 8 



Roots and relative normalization constants C^/C^ for the 10 B nucleus. 





(4a 1/2 )(4p 1/2 )(2p3 /2 ) 


(4s 1/2 )(6p 3/2 ) 




± V v J Ulo 


CM) /(7(°) 


T? nnt's 

J. IA ' ULo 


(7( £ ) /(7(°) 


VV 1 V T 1 \ 1 j 1 \, 


-2.4267116231 
-1 2816012954 


1.0 


-2.4267116231 
-1 2816012954 


1.0 


BARTLET 

1 J J 1111 UJ J -L 


-1.0057745600 
1.5278733354 


-0.6712940280 


-1.0066449150 
1.6322400952 


-0.5602373564 


HEISENBERG 


-0.9926668226 
0.5786209032 


-0.8464804834 


-0.9926668226 
0.5786209032 


-0.8464804834 


MAJORANA 


-0.8743416859 
4.6270236827 


-0.5924519571 


-0.8590657779 
0.2787982623 


-2.1204747128 
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Table 9 



Binding energies and rms radii for the B nucleus 



Potential 


(4s 1/2 )(4p 1/2 ) (2p 3/2 ) 


(4si/ 2 )(6p 3 / 2 ) 




EAA UAA 


EAA UAA 


Bl[29] 


54.787 54.542 
(2.699) (2.703) 


41.955 41.615 
(2.753) (2.758) 


MS3[30,31] 


33.475 32.941 
(2.648) (2.665) 


22.504 22.252 
( 2.756) (2.764) 


GPDT[32] 


17.928 17.827 
(3.139) (3.144) 


9.350 9.269 
(3.361) (3.369) 
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Table 10 

The J (I m , n) coefficient 



n 


1 


2 


3 


4 


5 


1 


16 


100 


400 


1225 


3136 


1 


-8 


-100 


-600 


-2450 


-7840 


2 


1 


35 


345 


1960 


8036 


3 




-5 


-290/3 


-2450/3 


-13328/3 


4 




1/4 


14 


1169/6 


4424/3 


5 






-1 


-329/12 


-4634/15 


6 






1/36 


161/72 


749/18 


7 








-7/72 


-161/45 


8 








1/576 


17/90 


9 










-1/180 


10 










1/14400 
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Table 11 

The I a p{^m,n) coefficient. The symbol [— n] denotes the exponential 10 



n 





1 


2 


3 


4 


5 


6 


7 





1 


4 


10 


20 


35 


56 


84 


120 


1 


























2 




1 


5 


15 


35 


70 


126 


210 


3 






-1 


-6 


-21 


-56 


-126 


-252 


4 






1/4 


7/3 


133/12 


112/3 


203/2 


238 


5 








-1/3 


-35/12 


-14 


-49 


-140 


6 








1/36 


35/72 


623/180 


15.9833 


56.583 


7 










-1/24 


-8/15 


-3.4833 


-15.833 


8 










1/576 


5.1389[-2] 


0.51458 


3.1268 


9 












-2.777[-3] 


-5.0463 [-2] 


-0.43638 


10 












6.9444[-5] 


3.1713[-3] 


4.2774[-2] 


11 














-1.1574[-4] 


-2.8770[-3] 


12 














1.9290 [-6] 


1.2676 [-4] 


13 
















-3.3069[-6] 


14 
















3.9367[-8] 



72 



